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Abstract. We characterize asymptotic collective behavior of rectangular random matrices, 
the sizes of which tend to infinity at different rates. It appears that one can compute the 
limits of all non commutative moments (thus all spectral properties) of the random matrices we 
consider because, when embedded in a space of larger square matrices, independent rectangular 
random matrices are asymptotically free with amalgamation over a subalgebra. Therefore, we 
can define a "rectangular free convolution", which allows to deduce the singular values of the 
sum of two large independent rectangular random matrices from the individual singular values. 
This convolution is linearized by cumulants and by an analytic integral transform, that we called 
the "rectangular J?-transform" . 

Introduction 

The first problem we are going to deal with in this paper is the modeling of asymptotic 
collective behavior of independent rectangular random matrices. In order to explain to the 
reader the way we will treat this problem, let us recall him the work already done by Wigner, 
Pastur, Marchenko, Girko, Bai, Voiculescu,... for square random matrices. First of all, in the 
50's Wigner considered self-adjoint random matrices with gaussian entries (GUE) and proved 
that the spectral law (i.e. the uniform distribution on the set of eigenvalues) of these random 
matrices converges to the so-called semicircle law. This result was improved, and other results 
giving the asymptotic spectral law of random matrices were proved (see, among many other 
sources, |PL02| ). In the same time, people studied the local structure of the spectrum of random 
matrices (see, e.g. |M67] ). but this is not the kind of problem we are going to study here. A new 
point of view was adopted in the early 90's by Voiculescu, who proposed a way to compute the 
asymptotic normalized trace (when the dimension n of the matrices goes to infinity) of products 

M{sunY'---M{sk,ny^ 

of random matrices taken among a family M{l,n), M{2,n), . . . of independent random nxn 
matrices and their adjoints, the only hypothesis being a bound on the norms of the matrices, 
the fact that the matrices M{l,n), M(2,n), . . . have all a limit singular law (the singular law 
of a rectangular matrix is the uniform distribution on the set of its singular values, i.e. of the 
eigenvalues of the absolute value of the matrix) , and a property of invariance of the distributions 
under an action of the unitary group. These results can be deduced from the fundamental article 
[V91j . but are presented under this form in [HPOOj . The advantage of being able to compute 
the limit of such normalized traces is that it gives us the asymptotic normalized trace of any 
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noncommutative polynomial in our independent random matrices. Hence, since the normalized 
trace of the fc-th power of a matrix is the A:-th moment of its spectral law, we are able to give 
the asymptotic singular law of any polynomial of our random matrices. This is why this work 
is said to model the asymptotic collective behavior of independent square random matrices. For 
example, it can be proved (combine results of [HPOOj and |HL00j ) that the asymptotic singular 
law of the sum of two independent random matrices whose distributions are invariant under left 
and right actions of the unitary group and whose singular laws converge weakly to probability 
measures ^i, fi2, only depends on fj.i and fi2, and can be expressed easily from fii and fi2'- it 
is the probability measure on [0, oo), the symmetrization of which is the free convolution of the 
symmetrizations of /xi, fi2- We often present the similar result for hermitian matrices, for which 
we work with spectral law in the place of singular law, and for which no symmetrization is 
necessary, but here, we shall work with rectangular matrices, which cannot be hermitian, so the 
square analogue of our work will be found in non hermitian matrices. In this text, as in the work 
of Voiculescu presented above, we will propose a way to compute the asymptotic normalized 
trace of products of random matrices taken among a family of independent rectangular random 
matrices, whose sizes tend to infinity, but with different rates. The notion involved, similarly 
to freeness in Voiculescu's modeling of asymptotics of square random matrices, is freeness with 
amalgamation over a finite dimensional subalgebra. This notion arises from operator-valued 
free probability theory, but we chose to use the point of view of operator-valued free probability 
only when necessary, because this point of view is not satisfying in all cases: if one uses this 
point of view, then he has to consider the case where the ratio of sizes of our random matrices 
tends to zero separately. Moreover, the point of view of operator-valued free probability in the 
asymptotics of rectangular random matrices is developed in another paper, where we analyze 
the related free entropy and Fischer information ( |B-G2| ). 

This modeling of asymptotics of rectangular random matrices will allow us to define, for 
A G [0,1], a binary operation on the set of symmetric probability measures, called free 
convolution with ratio A, and denoted by m^. For /Ui, /X2 symmetric probability measures, fiim^fi2 
is defined to be the limit of the singular law of a sum of two independent rectangular random 
matrices, whose dimensions tend to infinity in a ratio A, one of them being bi-unitarily invariant, 
and whose singular laws tend to /ii,/X2. In the second part of this paper, after having analyzed 
related cumulants, we construct an analytic integral transform which linearizes Ea;^ (like the R- 
transform does for free convolution) and give examples. This part of this paper is the base of 
other researches. Firstly, in [B-Gl], we study the related infinite divisibility: it is proved that the 
set of ffl;^-infinitely divisible distributions is in a deep correspondence with the set of symmetric 
classical infinitely divisible distributions. Secondly, in |BBG07] . we study some questions related 
to the support and the regularity of Him^fj,2- 

Acknowledgments. We would like to thank Philippe Biane, Dan Voiculescu, and Piotr 
Sniady for useful discussions, as well as Thierry Cabanal-Duvillard, who organized the workshop 
"Journee Probabilites Libres" at MAP5 in June 2004, where the author had the opportunity 
to have some of these discussions. Also, we would like to thank Cecile Martineau for her 
contribution to the english version of this paper. 

1. Asymptotic behavior of rectangular random matrices 

1.1. Notes on reduction of rectangular random matrices. In this subsection, we recall 
some simple facts about polar decomposition of rectangular complex matrices, which can be 
found in [HJ91] . Consider a pxq matrix M. When p < q, we will denote by \M\ the only pxp 
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positive hermitian matrix H such that we can write M = HT, with T a pxq matrix such that 
TT* = Ip. In this case, the matrix \M\ is the square root of MM* . When p > q, we will denote 
by \M\ the only qxq positive hermitian matrix H such that we can write M = TH, with T a pxq 
matrix such that T*T = Iq. In this case, the matrix \M\ is the square root of M*M. In both 
cases, the spectrum of \M\ is the only (up to a permutation) family (hi, /12 • • •) of nonnegative 
real numbers such that one can write 

M = U[6ihi]i<i<pV, 

l<j<q 

with U, V respectively pxp, qxq unitary matrices. These numbers are called singular values 
of M. The uniform distribution on /ii, /12... will be called the singular law of M (whereas the 
uniform distribution on the spectrum of an hermitian matrix is called its spectral law). Note 
that in both cases, for any a > 0, 

Tr \M\" = Tx{MM*Y^ = Tr(M*M)"/^ 



mm{p, q) min(p, q) min(p, q) 

is the a-th moment of the singular law of M. 

A random matrix is said to be hi-unitarily invariant if its distribution is invariant under the 
left and right actions of the unitary group. 

1.2. Definitions. For definitions of singular values, singular law, spectral law, (bi-)unitarily 
invariant and uniform random unitary matrices, we refer to subsection II. li For all positive 
integer d, we denote the set {l,...,d} by [d]. We denote the normalized trace of a square 
matrix X by tr X (whatever the size of X is) . 

Consider a positive integer d, and sequences qi{n) , . . . , qdin) of pairwise distinct positive 
integers which all tend to infinity as n tends to infinity, and such that for all k € [d], 

(1) > Pk>0, 

n n^oo 

at most one of the p^'s being zero. We are going to work with rectangular random matrices 
of sizes qk (n) x qi (n) . Free probability's modeling of the asymptotic behavior of square random 
matrices relies on a comparison between random matrices and elements of an algebra arising 
from operator algebra theory. To product analogous results for rectangular random matrices, we 
have embedded our matrices of sizes qk{n)xqi{n) {k,l S [d]) in an algebra. Let us assume that 
for all n, qi{n) + • • • + qd{n) = n (since what we will prove would obviously also work when n is 
replaced by a subsequence, it is not a real restriction). Now, nxn matrices will be represented 
as d'^ block matrices, such that for all k,l £ [d], the {k,l)-th. block is a qk{n)xqi{n) matrix. 

For all k, I G [d\, for all qk[n)xqi{n) matrix M, let us denote by M the "nxn extension of M", 
that is the nxn matrix with (i,j)-th block M if (i,j) = {k,l), and zero in the other case. Note 
we have preservation of adjoints (M* = M*), and of products: 



\/k,l,k',l' £ [d],yM of size qk{n)xqi{n),\/N of size qk>{n)xqii{n),MN 

At last, let us define, for all k £ [d\, the projector 

Pk{n) ■■= Iqk{n)- 




if / = k', 

in the other case. 



Let us formalize the structure we inherit. Consider a *-algebra A endowed with a family 
{pi, . . . ,pd) of non zero self-adjoint projectors (i.e. 'ii,pf = Pi) which are pairwise orthogonal 
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(i.e. yi / j,PiPj = 0), and such that pi + • • • + Pd = 1- Any element x oi A can then be 
represented 

xn ■■■ xid 
x= \ ': , 

_Xdl ■ ■ ■ Xdd_ 

where \/i,j^Xij = Pixpj. This notation is compatible with the product and the involution. 
Assume each subalgebra PkApk to be endowed with a tracial state (i.e. fkiPk) = 1 and for 
all x,y £ pkApk, ^k{xy — yx) = 0) such that for all A;, Z € [d], x G pkApi, y £ piApk, 

(2) pk(pk{xy) = pmiyx), 

where {pi, . . . , pd) is still the sequence of nonnegative real numbers defined by ([1]). Recall that 
at most one of the pkS is zero. 

Definition 1.1. Such a family {A,pi, . . . ,Pd,'^i, ■ ■ ■ ,^d) will he called a (pi, . . . , p^^) -rectangular 
probability space. Elements of the union of the PkApi 's (k, I £ [d]) will be called simple elements. 

Example 1.2. A is the *-algebra of nxn complex matrices, Pi,---,Pd o,re pi{n) , . . . , pd^n) 
previously defined, and each cpk is ^^7^ Tr. Then (0j is satisfied when each pk is replaced by 
qk{n)/n. 

Consider a (pi, . . . , p^) -rectangular probability space {A,pi, . . . ,pd, (pi, ■ ■ ■ , fd)- Denote by V 
the linear span of the pkS, then P is an algebra and the function E : ^ — > C^, which maps 
X £ A to E(x) = Ylk=i^k{xkk)Pk is a conditional expectation from A to V, i.e. E(l) = 1 and 
y{d,a,d') £V xAx V,E{dad') = dE{a)d' . 

The following definition, first appeared in |V95| , gives the right notion to describe asymptotics 
of independent random matrices. 

Definition 1.3. In an *-algebra B endowed with a *-subalgebra C and a conditional expectation 
Ec from B to C, a family {C C Ba)a&A of *-subalgebras is said to be free with amalgamation 
over C if for all m > 1, for all ai ^ ■ ■ ■ ^ am £ A, for all xi, . . . , Xm elements of respectively 
Bai , ■ ■ ■ , Bam one has 

(3) Ec(xi) = ••• = Ec(x™) = =^ Ec{xi---xm) = 0. 

A family of subsets of B is said to be free with amalgamation over C if the subsets are contained 
in *-subalgebras which are free with amalgamation over C. 

Definition 1.4. In a (pi, . . . , pd) -rectangular probability space, the P-distribution of a family 
{aj)j(zj of simple elements is the function which maps any polynomial P in the noncommutative 
variables {Xj,X*)j^j , with coefficients in V, to E{P{aj,a*)j(^j). 

The last notion we have to introduce is the convergence in D- distribution. Note that the P's 
of different spaces can all be identified with C^. 

Definition 1.5. • If for all n, {An,Pi,ni ■ ■ ■ ,Pd,n,^i,n, ■ ■ ■ ,'Pd,n) is a (pi,„, . . . , pd^n) -probability 
space such that 

iPl,n, ■ ■ ■ , Pd,n) > iPl,---,Pd), 

n— >oo 

a family {aj{n))j(zj of simple elements of An is said to converge in D-distribution, when n goes 

to infinity, to a family {aj)j(^j of elements of A if the T)- distributions converge pointwise. 

• In the context of example if for all n, aj{n)'s are nxn random matrices, convergence 
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in P-distribution in probability of the family {aj{n))j^j to {aj)j(zj is the convergence in prob- 
ability, when n oo, of the random variable E(P{aj{n),aj{n)*)j^j) to E(P(aj, a^)jgj) for all 
polynomial P in the noncommutative variables (Xj , X*)j^j . 

1.3. Statement of the theorems about random matrices. Let, for s G N, fc, / G [d], n > 1, 
R{s,k,l,n) be a qk{n)xqi{n) random matrix such that for all n, the family (R{s,k,l,n))s^k,l is 
independent. We suppose that for all s GN,k ^ I £ [d], R{s, k, I, n) is bi-unitarily invariant and 
that its singular law converge in probability to a deterministic probability measure. We also 
suppose that for all s € N odd, for all fee [d], R{s,k,k,n) is unitary uniform and that for all 
s € N even, for all k £ [d], R{s,k,k,n) is hermitian, unitarily invariant and that its spectral 
law converge in probability to a deterministic probability measure. Consider also, for n > 1, a 
family (C(i,n))jg7 of diagonal deterministic matrices, each C{i,n) having the size qki{n)xqi-{n), 
for a certain (fej, /j) G [d]x [d], which contains pi{n), . . . ,Pdin) and is stable under product and 
taking adjoints. We suppose that for all i £ I such that ki = li, the normalized trace of C{i,n) 
has a finite limit when n goes to infinity. 

We suppose that for the operator norm || • || associated to canonical hermitian norms these 
matrices are uniformly bounded (in n, s, k, I, i). 

The hypothesis of our main theorems are not the same if one supposes all /o^'s to be positive 
or not. In the following theorem, all p^s are supposed to be positive. 

Theorem 1.6 (Case where all p^'s are positive). Under these hypothesis, the family {R{s, k, I, n) ; 
s G N, A;, / G [c^]) U {C{i, n) ; i £ I) converges in D- distribution in probability to a family which is 
free with amalgamation over D. 

In the following theorem, we suppose one of the p^'s, say pi, to be zero. We supposed that 
at most one of the p^s is zero, so P2, ■ ■ ■ , Pd are positive. 

Theorem 1.7 (Case where pi = 0, other p^'s are positive). Assume moreover that for all s,k ^ 
/, n, R{s, k, /, n) has deterministic singular values and that for all s £ 'N even, k £ [d], R{s, k, k, n) 
has a deterministic spectrum. Then the family (R{s,k,l,n) ; s £ f^,k,l £ [d]) U {C{i,n) ; i £ I) 
converges in D- distribution in probability to a family which is free with amalgamation over T>. 

Corollary 1.8 (Product of rectangular matrices with asymptotic null ratio). Let qi(n),q2{n) 
be sequences which tend to infinity in such a way that qi{n) / q2{n) tends to zero. If, for all n, 
Mi{n),M2{n) are independent qi{n)xq2{n) bi-unitarily invariant rectangular random matrices 
with deterministic singular values and limit singular laws, then for all p > ii, ■ ■ ■ ,ip £ {1, 2}, 

• Mjj (n)*Mj2(n) • • • Mi^^_^{n)* Mi,^^{n) has null limit normalized trace, 

• Mi-^(n)* Mi^(n) ■ ■ ■ Mi^^_j^{n)* Mi2^{n) has null limit normalized trace, except if any matrix 
Mi{n) or M2{n) in the product is followed by respectively Mi{n)* or M2{n)* , in which case the 
limit can be computed using the well known asymptotic freeness of square random matrices. 

1.4. Comparison with already existing results. There already exists some results which 
allow to compute asymptotics of normalized traces of products of rectangular random matrices: 
on one hand the results about Wishart matrices and on the other hand the theorem 4.1 of |Sh96] . 

Theorems 11.61 and 11.71 allow us to compute the limit (for convergence in probability) of all 
normalized traces of matrices of the type 

(4) R{si, ki,li,nY^C{ii,n)R{s2, k2,l2,nf^C{i2,n) ■ ■ ■ C{ip^i,n)R{sp, kp, Ip, nfpC{ip, n), 
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with p > 1, Si, . . . ,Sp G N, ki,li, . . . ,kp,lp £ [d], ei,. . . ,ep £ {1, *}, and h,. ■ ■ , ip-i G / such 
that the product is possible and is a square matrix. 

The Wishart matrices are matrices of the type RR* , where i? is a rectangular gaussian matrix. 
In |HP00j , |CC04j or |CDM05] , the authors give results of asymptotic freeness for such random 
matrices. In the case we are interested in, where the random matrix R is bi-unitarily invariant, 
this asymptotic freeness is only a consequence of the asymptotic freeness of unitarily invariant 
hermitian random matrices. With our notations, it does not allow to suppose that pi = 0, and 
to consider products where rectangular matrices are not immediately followed by their 
adjoint, as the one of 

Theorem 4.1 of the paper [Sh96] by Shlyakhtenko allows, embedding rectangular matrices 

R 

to consider products of the type of the one of 



in larger self-adjoint matrices {R 



R* 



([¥]), but only when the pj's are positive, the matrices R{s,k,l,n) are gaussian, and the 
deterministic matrices are diagonal and have, in a certain sense, a limit in L°°[0, 1]. Moreover, 
this theorem states the convergence of the expectations of the normalized traces, but not their 
convergence in probability (even though this could it is possible to expect to pass, using the 
concentration results of i( iZOOj . from convergence in expectation to convergence in probability 
in Shlyakhtenko's result). 

As an example, no general result had been proved about the convergence, when n,p — > oo such 
that ^ ^ A > 0, of the singular law of M {n, p) + N (n, p) , where M(n,p),N{n,p) are independent 
bi-unitarily invariant nxp random matrices with limit singular laws. We did not even know if this 
limit existed, even though some computations have been done in some particular cases ( |DS04] . 
|HLN04j . |HLN05j )). Indeed, the computation, by the moments method, of the singular law of 
M(n,p) + N{n,p) involves the computation of the normalized trace of 

[{M{n,p) + Nin,p)){M{n,p) + iV(n,p))*]^ 

the expansion of which contains products where some rectangular matrices are not followed by 
their adjoints. 



1.5. Proof of theorem 11.61 (case where all p^'s are positive). Let us introduce, on the 
space of nxn matrices, the Schatten p-norms (p G [l,oo]), defined, for p < oo, by 

\\X\\p = {tT\X\P)K 

and 1 1. 1 loo = I Ml- They satisfy classical Holder inequalities ( [N74| ) . 

Step I. Fix k,l £ [d] and consider, for n > 1, a bi-unitarily invariant qk{n)xqi{n) random 
matrix R{n) such that the sequence {R{n)) is uniformly bounded for ||.|| by C > and such that 
the singular law of R{n) converge in probability to a deterministic probability measure. Note 
first that, by Carleman criterion, it implies that its singular law converges in probability to a 
probability measure /i. It is easy, transferring the proof of lemma 4.3.10 p. 160 of [HPOO], to 
prove that the distribution of R{n) can be realized as the distribution of U{n)A.{n)V{n) where 
U{n), A(n), V{n) are independent, U{n), V{n) are respectively qk{n)xqk(n), qi(n)xqi(n) uniform 
random unitary matrices, and A(n) is a rectangular diagonal matrix, 

A(n) = [SlX{i,n)]i<i<q^(^n) with < A(l,n) < • • • < X{mm{qkin), qi{n)),n). 

l<l<qi{n) 
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Choose deterministic < .^(l,n) < ••• < S,{mm{qk{n),qi{n)),n) < C such that the uniform 
distribution on the ^(l,n)'s converges weakly to /i. Now set 

E{n) := [(^■A(i,n)]i<j<q^„) and N{n) := U {n)E{n)V (n) . 
i<i<g;(") 

Let us prove that for ah P polynomial in noncommutative random variables X,X*, for all 
p G [1,cxd), \\P(R(n)) — P{N{n))\\p converges in probability to zero. 

By lemma 4.3.4 p. 152 of [HPOOj . for all p>l, 

\\R{nr-Ninr\\, = ||i?(n) - iV(n)||, = ||A(n) - 

converges in probability to zero. So, for all m > 1, ei, . . . , £ {-j *}j 

I |^(n)^i • • • R{nY^ - N{nY^ ■ ■ ■ N{nY^ \ \p 

< YlZi ■ ■ ■ RinY^-^{R{nY^ - N{nY^)N{nY'+^ ■ ■ ■ iV(n)^'"||p 

< ET=iC^-'mnY'-N{nY')\\p, 

which proves the result. 

Note that in the case where i? is a uniformly bounded unitarily invariant qki'n)xqii.{n) hermitian 
random matrix, the moments of the spectral law of which converge in probability, then the same 
work can be done, replacing V{n) by U{nY and singular values by eigenvalues, and the same 
conclusion holds: for all P G C[X], for all p G [l,oo), \\P{R{n)) — P{N{n))\\p converges in 
probability to zero. 

Step II. By hypothesis, there exists ii, . . . , G / such that for all n, k, C{ik,n) = Pk{n) and 
there exists a map (3 from I^ to / and an involution * from / such that for all n > 0,i,j G / 
such that li = kj, 

(5) Cii,n)Cij,n) = C^, j),n), Cii,nY = C(*(i),n). 

Consider, in a (pi, . . . , -rectangular probability space {A,pi, . . . ,pd, fi, ■ ■ ■ , (fd), a family 
a{s,k,l) {s G f^,k,l G [d]), c{i) (i G I) of elements of A, with c{ii) = pi,.-., c{id) = Pd, 
whose D-distribution is defined by the following rules: 

(i) For all s G N, fc, / G [d], a{s, k, I) G pkApi, and for all i £ I, c{i) G PkiApi^. 

(ii) For all i,j G / such that k = kj, c{i)c{j) = c{P{i, j)), c{iY = c(*(i)). 

(iii) For all s £ n,k ^ I £ [d], for ah r > 0, 

ifk[ici'{s,k,l)a{s,k,iyY']= lim tr{R{s,k,l,n)R{s,k,l,nYY , 

n—>oo 

ifi[{a{s,k,lYa{s,k,lY^] = li™ tv{R{s, k,l,nY R{s, k,l,n)Y . 

n— »oo 

(iv) For all s G N odd, k G [d], a{s,k,k) is unitary in p^Apk, and for all r G Z — {0}, 
iPk{a{s, k, kY) = 0. 

(v) For all s G N even. A; G [d], a(s, k, k) is self-adjoint and for all r > 0, 

(pk[a{s, k, kY] = lim tr R{s, k, k,nY , 

n^oo 

(vi) For alH G / such that ki = li, ipk^^c^iY is the limit of the normalized trace of C{i,n). 

(vii) The family {a[s, k, ?)} (s G N, fc, / G [d]), {c(i) ; z G /} is free with amalgamation over T). 

Note that such a space and such a family exist, they are given by the free product with 
amalgamation r |VDN91j . [S98] ). 
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Let US denote by C the linear span of the c(z)'s. By ([5]) and (ii), for all n, there is a *-algebra 
morphism from C to the linear span C„ of the C{i,nys. By (vi), for all c G C, for all k G [d], 
the normalized trace of the k-th diagonal block of ipnic) tends to (pk{ckk)- 

Since 1 G C, it suffices to prove that for all m > 0, for all {sj,kj,lj)j^^^^ G (Nx [djx [d])'", for 
all Pi, . . . ,Pm polynomials in the noncommutative variables X, X*, for all (i(0), . . . ,d{m) G C, 
EinjLi PjWsj, kj,lj,n))'ipn{d{j))] converges in probability to E[d(0) YYjLi Pj{aisj,kj,lj))d{j)]. 

Step III. By the formula Tr XY = TrYX and since all p^s are positive, it is easy to see that 
we can remove d{0). 

Step IV. Then, we claim that it suffices to prove it when the R{s, k, I, n)'s are replaced by the 
N{s,k,l,nys, where the A^(s, /c, /, n)'s, are defined from the /c, /, n)'s like N{n) from R{n) 
in the step I. Indeed, the normalized trace of any diagonal block of 

m m 

J] Pj{R{sj, kj,lj,n))MdU)) - n Pji^i^J^ kj,lj,n))Md{j)) 

has an absolute value less or equal than the || • ||i norm of their difference times one of the 
n/qk{n), which is less or equal, by a decomposition like in step I, than 

m 

^^^^y,C'^D'^-'\\Pj{Ris,,k,J,,n)) - P,{N{s,,k,,l,,n))\\,, 
fce[c(] qk{n) ^ 

where D = max \Pj{z,z')\ which converges in probability to zero, by step I and because all 

je[m] 

PkS are positive. 

So let us prove that for all k £ [d], the normalized trace of the k-th diagonal block of 

m 

l[Pj{N{s,,k,,l,,n))Md{j)) 
i=i 

converges in probability to the k-th coordinate of E[]^^j^ Pj{a{.Sj, kj, lj))d{j)]. 

Step V. We will prove it by induction on m. If m = 0, the result is clear. In the other case, let 
us suppose the result to be proved to the ranks 0, . . . , m— 1. Let us denote x{j) = Pj{a{sj, kj,lj)) 
and X{j,n) = Pj{N{sj,kj,lj,n)). Since C is an algebra containing pi, . . . ,P(i, if, for a certain 
j G [m], one would replace X{j, n) by one of the pk{n)'s and x{j) by the corresponding p^, then 
the result would follow from the induction hypothesis (and from step III if j = 1). Thus, by 
linearity, one can, for all j G [m], add a linear combination of the pfc(n)'s to X{j,n) and the 
same linear combination of the p^s to x{j). Therefore we can assume that for all j, E{x{j)) = 0. 
By linearity, we can moreover suppose that for all j, x{j) is a simple element. Note that these 
suppositions imply that for all n,j, X(n,j) has only one non zero block, E(X(j, n)) converges 
in probability to zero. 

From now on, we use the cyclic order on [m]. This means that to put the index m + 1 on an 
element amounts to put the index 1. As simple elements, the x(j)'s have the following property: 
the product of any of the x(j)'s by one of the p^s is either zero or x{j), and the same holds for 
n)'s with pfc(n)'s. Thus, if, for a certain j G [m] such that {sj,kj,lj) = (sj+i, fcj+i, Zj+i), 
one would replace d{j) by a linear combination of the pfc(n)'s, then the result would follow from 
the induction hypothesis (and from step III if j = m). So one can suppose that for all j G [m — 1] 
such that {sj,kj,lj) = (sj+i, fcj+i, /j+i), E{d{j)) = 0. 
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At last, for j £ [m] such that (sj, kj, Ij) ^ {sj+i, fcj+i, ^j+i), one can write 
d{j) = d!{j) + \k{j)Pk, with E(d'(j)) = 0. 

ke[d] 

Then, by hnearity, one can suppose that for ah such j, E{d{j)) = or 3/;: € [d], d{j) = p^. 

To conclude, we have shown that we only have to prove that E[J3J1-^ n)-(/'n(c^(j))] con- 
verges in probability to E[]^™^]^ under the hypothesis that all x(j)'s are simple elements 
and satisfy E(2;(j)) = 0, and for all j G [m], 

E{d{j)) = or [d{j) G and {sj,kj,lj) 7^ (s^+i, /j+i)]. 

Note that under this hypothesis, by freeness with amalgamation over D, E[]^™^]^ = 0. 

Step VI. For all s, k, I, n with k ^ I, the matrix N{s, k, I, n) has been defined by the work of 
step I, so we can write 

N{s, k, I, n) = [/(s. A;, /, n)H(s, fc, /, n)y(s, fc, /, n), 

with [/(s, /c, /, n), V{s, k, I, n) uniform random unitary matrices with respective sizes Qkin), qi{n), 
and 3(3, k, I, n) a deterministic diagonal qk{n)xqi{n) matrix. In the same way, when s is even, 
N{s, k, k, n) can be written 

A^(s, A;, l,n) = U{s, k, k, n)H(s, k, k, n)U* (s, k, k, n), 

and when s is odd, N{s, k, k, n) is a uniform random matrix of Ufc(n). 

Consider j G [m]. If kj 7^ /j, then X{j, n) is of one of the 4 following types: 

(1) U{sj,kj,lj,n)AjV{sj,kj,lj,n), where Aj is a {n)xqi.{n) diagonal matrix whose di- 
agonal terms are the images of the diagonal terms of H(sj, kj, lj,n.) by an odd polynomial. 

(2) V{sj,kj,lj,n)*AjU{sj,kj,lj,n)*, where Aj is a qi.{n)xqk.{n) diagonal matrix whose di- 
agonal terms are the images of the diagonal terms of H(sj, kj, lj,n) by an odd polynomial. 

(3) U{sj,kj,lj,n)AjU{sj,kj,lj,n)* , where Aj is a qkj{'n)xqk^{n) diagonal matrix whose di- 
agonal terms are the images of the diagonal terms of H(sj , kj , Ij , n) by an even polynomial. 

(4) V{sj, kj,lj,n)*AjV{sj, kj,lj,n), where Aj is a g/ - {n)xqi . (n) diagonal matrix whose diag- 
onal terms are the images of the diagonal terms of H(sj, kj,lj,n) by an even polynomial. 

And if kj = Ij, when Sj is even, only the last type is possible, and when Sj is odd, X{j,n) is 
a linear combination of positive or negative powers of N{sj,kj,lj,n). By linearity again, this 
combination can be replaced by a non zero power of N{sj, kj,lj,n). 

Note that in all types previously presented, the constant diagonal matrix A,- satisfies E(Ao) — 

0. Moreover, for all n, the family of random unitary matrices arising in these decompositions 
are independent, so the result follows from case 1°) of proposition 14.41 and remark [4.51 

1.6. Proof of theorem 11.71 (case where pi = 0, all other p^'s are positive). We shall 
follow closely the proof of theorem 11.61 (case where all p^s are positive) given in the previous 
section, but because of the weaker hypothesis (singular values of non hermitian random matrices 
and spectrum of hermitian ones are now supposed to be deterministic), we shall skip steps I 
and IV, in which non hermitian (resp. hermitian) random matrices where approximated by 



10 



FLORENT BENAYCH-GEORGES 



random matrices with deterministic singular values (resp. deterministic spectrum). However, 
the hypothesis pi = will make things slightly harder. 

Copying step II of the previous proof, we have to prove, for all m > 0, the following propo- 
sition: H{m):= "for all {sj,kj,lj)j^[fn] £ (Nx [djx [d])"^, for all Pi,...,Pm polynomials in the 
noncommutative variables X,X* , for all d{0), . . . , d{m) G C, for all A; G [d], 

m 

E[MdmllPjiR{sj,kj,lj,n))i^nidim 
i=i 

converges in probability to E[(i(0) YlJLi ^ji^i^j^ 

Note first that by linearity, it suffices to prove it when d{0),d{m) are supposed to be simple, 
which implies that there exists u,v £ [d] such that Pud{0) = d{0), d{m)pv = d{m). Let us denote 
by H{m,u,v) the proposition H in the particular case where Pud{0) = d{0), d{m)py = d{m). 
Let us also denote by H'[m) the proposition H[m) in the particular case where (i(0) = 1. 

Note first that H{m, u, v) is obvious when u ^ v. Moreover, like in step III of the previous 
section, we see that for all u £ [d] — {1}, H'(m) implies H{m,u,u). Thus we only have to prove 
that for all m > 0, H'{m) and H{m, 1, 1) hold. 

Let us prove it by induction on m. For m = 0, the result is immediate. Consider m > 1 such 
that H'{1) and H{1, 1, 1) hold for all I < m. 

Let us first prove H'{m). Consider (sj, kj, lj)j^[rn] ^ (^^^ M>< [^])™) Pi-, ■ ■ ■ ■, Pm polynomials in 
the noncommutative variables X,X* , . . . ,d{m) £ C. Let us prove that 

m 

E[l[Pj{R{sj,k,,lj,n))Md{M 

converges in probability to E[]^J1-^ Pj(a(sj, A;j, Let us denote x{j) = Pj{a{sj,kj,lj)) 

and X{j,n) = Pj{N{sj,kj,lj,n)). Similarly to step V of the previous section, by the induction 
hypothesis and the fact that H'{m — 1) implies H{m — 1, u, u), we can suppose that for all n, j, 
X(n,j) has only one non zero block and E(X(j, n)) converges in probability to zero. We can 
also suppose that for all j G [m] such that with the cyclic order, {sj,kj,lj) = (sj+i, /cj+i, /j+i), 
E{d{j)) = or there exists k G [d\ such that d{j) = p^. 

Under these assumptions, we only have to prove that ^lYYjLi X (j , n)'ipn{d{j))] converges in 
probability to E[J3J1^ x{j)d{j)] = 0. Following step VI of the previous section, it appears to be 
an application of case 1°) of proposition 14.41 and remark [4.51 

Similarly, H{m, 1, 1) is an application of case 2°) of proposition 14.41 

2. The rectangular free additive convolution 

For I' probability measure on R, denote by P the symmetrization of z^, which is the probability 
measure defined by 9(B) = ^{u{B) + u{-B)) for all Borel set B. 

One of the interests of modeling the asymptotic behavior of square random matrices with free 
probabilities is the possibility to compute the asymptotic spectral law of an hermitian matrix 
which is a function of several independent random matrices. In particular ([V91]), if M(l,n), 
M(2, n) are independent hermitian unitarily invariant random matrices whose spectral laws 
tend respectively to ui , 1/2 when their dimension n goes to infinity, then the spectral law of 
M(l,n) + M(2,n) tends to z^ifflz^2 (free additive convolution of i'i,i'2), and the spectral law of 
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y^M(T~nyM(2, n)y^M(Tyn) (when the matrices are positive) tends to 1/1^1/2 (free multiphcative 
convolution of ^'1,^2)- In the same way (combine Theorem 4.3.11 of |HP00j and Propositions 
3.5, 3.6 of [HLOO] ). if the matrices M{l,n), M{2,n) are square, bi-unitarily invariant (and still 
independent), with asymptotic singular distributions I'l, V2, then the symmetrization of the 
singular distribution of M(l, n) + M(2, n) tends to ^101^2) and the push-forward, by the function 
X — > x^, of the singular distribution of M(l,n)M(2,n) is the free multiplicative convolution of 
the push-forwards, by x ^ x^, of z^i, 1/2- 

Now, suppose that M(l,n), M(2,n), instead of being square, are rectangular (independent 
and bi-unitarily invariant) qi{n)xq2{n) random matrices whose singular laws tend to z^i,z^2- We 
keep the notations introduced in section [L2l but suppose that d = 2: 

qi{n) Q2in-) 
qi{n) + q2{n) = n, qi{n),q2{n) — > 00, > Pi > 0, > p2 > 0. 

n— >oo n n^oo n n—>oo 

Since for 1 < q < q' , the singular law of a qxq' matrix M is related to the singular law u' of 
M* by z/' = {q' — q)/q'6o + q/q'i^, we can, without restriction, suppose that qi{n) < q2{n). 

First, the singular distribution of M(l, n)M(2, n)* tends to a distribution that can be com- 
puted by free probability theory. Indeed, for all r > 1 even, the r-th moment of the singular 
distribution of M{l,n)M{2,n)* is 

Tr(M(l,n)M(2,n)*M(2,n)M(l,n)*)'''^ _ g2(n) Tr(M(2, n)*M(2, n)M(l, n)*M(l, n))''/^ 
qi{n) qi{n)q2{n) 

which tends to P2I P\ (when p\ > 0) times the r/2-th moment of tiht2, where t\^T2 are the 
limit spectral distributions of M(2, n)*M(2, n), M(l, n)*M(l, n). t\^T2 can easily be computed 
from 1^1,1^2- Thus, if p\ > 0, the limit of the spectral distribution of M{l,n)M{2,n)* can be 
easily computed using tools of free probability. If pi = 0, then the explicit computation of the 
expectation and of the variance of ([6]) for r = 2 gives us the convergence of the singular law of 
M(l, n)M(2, n)* to the Dirac measure in zero. 

But until now, it wasn't possible to compute the singular law of M(l,n) -|- M(2,n) with the 
tools of free probability theory. Indeed, for r > 1 even, its r-th moment is given by 

^— Tr((M(l, n) + M(2, n))(M(l, n)* + M(l, n)*))"^, 

which can be expanded into the sum, over all i G {1,2}'', of 

(7) Tr M(ii , n)M(i2,n)* ■ ■ ■ M(v„i, n)M(v, n)*. 

Asymptotics of normalized traces like in ([7]) cannot be computed with free probability the- 
ory, because our random matrices are not square, and we cannot reduce the problem, as for 
M(l, n)M(2, n)*, to a problem which involves only independent square random matrices. Thus 
we have to use theorem 11.61 (or corollary 11.81 if pi = 0) to compute the asymptotic singular law 
of M(l,n) + M(2,n). 

Proposition-Definition 2.1. Let pi,pL2 be two compactly supported symmetric probability mea- 
sures on the real line. 

(a) Let, for all n > 1, M{l,n), M{2,n) be independent bi-unitarily invariant qi{n)xq2{n) ran- 
dom matrices, with deterministic singular values if pi = 0, uniformly hounded for \ \.\\, and such 
that for all i = 1,2, the symmetrization of the singular law of M{i, n) converges in probability to 
pLi. Then the symmetrization of the singular law of M{\,n) + M[2,n) converges in probability to 
a symmetric probability measure on the real line, denoted pia^fi2, which depends only on pi,p2, 
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and A := lim„_^+oo qi{n) /q2{n) ■ 

(h) /iifflA/^2 is the unique symmetric measure /i such that for all r €zN even, the r-th moment of 
fi is ipi{{{ai +a2)(ai + a2)*Y^), where ai, 02 are free with amalgamation over V elements of 
PiAp2, {^^Pi-,P2-,Vi-,'^2) is a (pi, P2) -probability space, and for all i = 1,2, for all r £ even, 
(pi{{aia*Y^) is the r-th moment of pi. 

(c) The support of fj,is^p2 is contained in the sum of the convex hulls of the supports of pi and 

^2- 

The binary operation on the set of compactly supported symmetric probabihty measures, 
called rectangular free convolution with ratio A, will be extended in section 13.91 to the set of 
symmetric probability measures, and the same result will stay true without any hypothesis of 
boundedness. 

Proof. By theorem 11.61 (or theorem 11.71 if pi = 0), if ai, 02 are elements of a {pi,p2)- 
probability space {A,pi,p2, ^1, ^2) as in (b), then (M(l, n), M(2, n)) converges in P-distribution 
in probability to (01,02). It implies that for all k > 1, the normalized trace of [(M(l,n) + 
M(2, n))(M(l, n) +M(2,n))*]'^ converges in probability to a finite and deterministic limit 
as n goes to infinity. 

Let us first prove that the m^'s are the moments of a probability measure with support on 
[0, (si + S2)^]. It suffices to prove that for any real polynomial P = c^X^ which is nonnegative 
on this interval, Ckm^ > 0. Note that since the m^'s only depend of the P-distribution of 
(oi, 02) one can suppose that for all n, for all i = 1,2, 

M{i,n) = U{i,n)[6iXi{k,n)]i<k<q^(ri)y{i,n), 

l<l<q2{n) 

where Xi{l,n), . . . ,\i{qi{n),n) are real numbers of the support of pi such that the uniform 
distribution on the Ai(A;,n)'s (/c = 1, . . . ,qi{n)) converges weakly to pi and U{i,n), V{i,n), are 
independent respectively qi{n)x qi{n), q2{n)xq2{n) uniform unitary random matrices. In this 
case, Ckm^ is the limit of the normalized trace of 

^ Cfe[(M(l, n) + M(2, n))(M(l, n) + M(2, n))*]^ 
k 

By hypothesis on P, this matrix is hermitian and nonnegative for all n, hence Ckm^ > 0. 

Since on the set of probability measures on a compact interval, convergence of moments is 
equivalent to weak convergence, (a) and (c) are proved. Moreover, (b) follows from the first 
paragraph of the proof. □ 

Remark 2.2. When A = 1, the rectangular free convolution with ratio A is the well known ad- 
ditive free convolution. To see it, apply the previous proposition to sequences of square matrices. 
We know (Theorem 4-3.11 of [HPOOj ) that independent bi-unitarily invariant square matrices 
are asymptotically free R-diagonal elements and that the symmetrization of the distribution of 
the absolute value of the sum of two free R-diagonal elements is the free convolution of the 
symmetrizations of the distributions of their absolute values (Proposition 3.5 0/ [HLOO] which 
proves that = m for A = 1. We will see later that the free convolution with ratio is also 
related to the free convolution. 
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3.1. Relation moments-cumulants in a (pi,/J2) probability space. For an introduction to 

the general theory of cumulants in a P-probability space and the particular case of a {pi , P2)) probability 
space, we refer to section 2 of |B-G2j . Basically, if {A,pi,p2,fi,<f2) is a (pi, /92)-probability 
space, then for all n, there are two linear functions Cn\cn^ from the n-th tensor product with 
amalgamation over D of „4 to C such that: 

(a) for all a G piAp2, for all n > 1, 

(8) Miaa*r) = Yl c[^^|(a® •••®a*) JJ cfj^|(a* • • • a), 

7reNC'(2n) VStt V&tt 

min V odd min V even 

(9) 9?2((a*ar) =^^11 c\^v\ia* ■ ■ ■ <^ a) JJ cf^^|(a • • • a*), 

7reNC'(2n) VStt V&tt 

min V odd min V even 

where NC'(2n) is the set of noncrossing partitions of [2n] where all blocks have an even cardi- 
nality. This formula can easily be deduced from equations (2.2), (2,6) and (b) of section 2.2 of 
FG2], 

(b) a G piAp2, for all n even, we have piCn\a a* ■ ■ ■ a*) = p2Cn\a* (8) a • • • a), 

(c) for all a,b piAp2 free with amalgamation over D, for all for all n even, we have 

(10) cl^\ia + 6) 8) (a + 6)* 8) • • • (a + b)*) = c^^\a ^ a* ■ ■ ■ ^ a*) + cl^\b ^ b* ^ ■ ■ ■ ^ b*). 

From this, we deduce, for A = pi/p2- 
Proposition 3.1. For a £ piAp2, n positive integer, 

(11) Miaa*r)= E A-^Wnciy^a), 

7reNC'(2n) V(^tt 

where for all r even, Cr{a) denotes ci^\a • • • a*), and for all n, e{n) denotes the number of 
blocks of IT with even minimum. 

3.2. The case A — 0, characterization of fflg- 

Lemma 3.2. If X = 0, then for all a E PiAp2, for all n > 1, C2n(o) is the n-th free cumulant of 
the element aa* of {piApi,^pi). 

Proof. It suffices to prove that for all n > 1, 

(12) Vi{{aa*T)= E n^2|y|(«)- 

7reNC(n) VStt 

By ([TT]) . the left hand term is equal to X^,reNC'{2n) e(7r)=o Ili/eTr "^IVI (^)- Let us introduce the set 
NCd(?^) of non-crossing partitions tt of [2n] such that for all k even in [2n], k — l^k. Note that 
NCd(2n) is contained in the set of elements of NC'(2n) in which no block has even minimum. We 
claim that the inverse inclusion is true. Indeed, suppose the existence of vr E NC'(2n)\ NCd(2n) 
in which no block has even minimum. Define k = min{A; E [2n] ; k even and k—l'^k}.k cannot 
be the minimum of its block in vr, and by minimality of k, its preceding element j in this block 
has to be even. Then, as vr is non-crossing, the set {j + 1,..., /c — 1} is a union of classes of tt, 
but its cardinality is odd, which is in contradiction with vr E NC'(2n). 
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Now, ()12p follows easily from the fact that NCd(2n) is in correspondence with NC(n) by the 
order-preserving bijection from NC(n) onto NCd(2n), that maps any noncrossing partition vr to 
the partitionthat links two elements i,j of [2n] if and only if the upper integer parts of i/2 and 
j/2 are linked by vr. □ 

Now we are able to give the link between the free convolution with null ratio ebq the free 
convolution a: 

Proposition 3.3. The free convolution with null ratio of two compactly supported symmetric 
probability measures is the unique symmetric probability whose push-forward by the square func- 
tion is the free convolution of the push-forwards by the square function of the two probabilities. 

Proof. Denote, for v probability measure, mn{i') its n-th moment, and ^(i^) its n-th free 
cumulant. Recall ( |S94j ) that free cumulants of probability measures are defined by the formula 

(13) Vn>l,m„(z.)= ^ H •^I^K'^)' 

7rGNC(n) BGtt 

and that the free convolution of vi , z/2 compactly supported probability measures is the only 
probability with free cumulants ^{vi) + ^n{^2) {n > !)• 

Consider two symmetric compactly supported probability measures fj,i, fJ-2- Denote by ui , 1/2 
their respective push- forwards by the square function. It suffices to prove that for all n > 1, 
"^2ri(A'iffloA*2) = 1 i-S- that for all n > 1, 

m2n(^lffloAt2) = ^ W[^\B\{^l) + ^\B\iy2))- 

7reNC{n) -BGtt 

Consider 01,02 G piAp2 free with amalgamation over T) such that 

= l,2,Vn > l,(/7i((aia*)") = m2n{l^i)- 

Then one has, for all n > 1, 

m2„(^iffloM2) = 92i[((ai + a2)(a* + 02))"] 

= n C|B|(ai + a2) bydHI) 

7reNCd(2n) -BGtt ^ ; ~^ , . 

But by lemma [221 for all i = 1,2, for all m even, Cm{ai) is the m/2-th free cumulant of the 
element a^a* of {Ai,(pi), i.e. the m/2-th free cumulant of Vi. So 

"l2n(/UlEBo^2) = ^ (•^IVl('^l) + •^|V|(Z^2)) = "ln(z^lEBl^2), 

o-eNC(n) veo- 

which proves the result. □ 

3.3. Generating series of the cumulants of an element of piAp2- We are going to derive 
from (jlip a formula which links generating series of the sequence (/9i((aa*)") and C2n(fl)- Note 
that such a formula could be derived from the very general Theorem 2.2.3 of |S98j . but the work 
needed to translate this result to our context and tu put it into the form we shall use it is as 
long as this section. 

In this section, we work in the field C((X)) of fractions of the ring C[[X]] of formal power 
series in X with complex coefhcients, endowed with the classical addition and product. Note 
that for any element S of C[[X]] with null constant term, F 1-^ F o S is a well defined operation 
on C{{X)). If moreover, the coefficient of X in is non null, then S has a left and right inverse 
for o, denoted by 5^^^^. 
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Lemma 3.4. Consider a sequence {c2n)neN* of complex numbers. Define the sequence {m^2l)nef^ 
by m^Q ^ = 1 and for each ra G N* , 

7reNC'(2n) Ben 

Define the formal power series 

C{X) = Y,C2nX\ m(^\x) = Y,mtlx\ 

n>l n>l 

Then we have C[X{\M^''\X) + 1){M^''\X) + 1)] = M^^\X). 

Proof. Step I. Define, for n positive integer and tt G NC'(2n), o(7r) to be the number of 
blocks of TT with odd minimum. Then we have 

7reNC'(2n) SGtt 

Indeed, if c denotes the cycle {2n —>-2n— >-2—>-l) of [2n], then c induces a permutation 
of NC'(2n) by 

: TT G NC'(2n) ^ R{^) := {c(B) ; 5 G 

For example, for vr = {{1, 4, 5, 6}, {2, 3}}, R{-k) = {{3, 4, 5, 6}, {1, 2}}. For tt G NC'(2n), for 
B block of TT, minc(5) is (minS) — 1 if minS > 1, and the second element of B minus 1 
in the other case. Note that in the first case, mini? and min c(i?) have different parities, 
whereas in the second one, they have the same parities, by definition of NC'(2n). Hence we have 
o{R{'k)) = e(7r) + 1. Thus 

A-s= E >^''^'^'\i^\B\= E A«^^^^^^ n ^i-i = E A°wn^i-i- 

7reNC'(2n) BEtt 7reNC'(2n) BStt 7reNC'(2n) Ben 

Step II. Thus we have to prove that C[X{A&\X) + 1){m'-''\x) + 1)] = M'^^\x), where 
= E„>i 45^", and = E ^^^^^ IT ^l^h 

7reNC'(2n) Ben 

Define, for each positive n and each j G [n], the "decomposition map" Dj, from the set of 
partitions of NC'(2n) in which the class of 1 has cardinality 2j to U NC'(2Zi) x • • • x 

h-\ \-l2j=n-j 

NC'(2/2j) (the set NC'(O) is considered as a singleton, on which e and a have value 0), which 
maps TT G NC'{n) to (tti, . . . ,7r2j) defined in the following way: if {1 = fci < . . . < k2j} is the 
class of 1 in tt, then for each r, let tTt be the restriction of tt to the interval \kr,kr+i[, with 
k2j+\ = 2n + 1. It is easy to see that the map Dj is well defined (recall that for all r, kr and 
kr+i have opposite parity), and that it is a bijection. Moreover, we have 

e(7r) = o(7ri) + e(7r2) H h o(7r2j_i) + e(7r2j), 

o(7r) = 1 + e(7ri) + o(7r2) H h e(7r2j-i) + o(7r2j). 

We denote by {X^)P the coefficient of X^ in a formal power series P in X. We have to show 
that for each n> 1, 

(X")C[X(m(°)(X) + 1)(m(^)(X) + 1)] = m£. 
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We have 



(X")C[X(M(°)(X) + 1)(M(^)(X) + 1)] 

n 
n 

5^c2,(x"-^')[(m(°)(x) + iy{M^'\x) + ly] 



(o) (e) (o) (e) 

Zli Zl2 Zl2j — 1 Zi2j 

j=l ZeN^i 

/iH \-l2j=n-j 

n 2j 

= E E A"(-H^(-)+-+^(-^) n n -1^1 

j=l ZgN^i 7rieNC'(2«i) r=lBe7Tr 

iiH \-l2j=n-j 

7r2ieNC'(2«2j) 

Hence the preHminaries about the bijections Dj tell us that 

(X")C[X(m(°)(X) + 1)(m('=)(X) + 1)] = ^ A<-) n ^1^1 = "^S- 

7reNC'{2n) -BStt 

□ 

Let us introduce two power series that will play a role in the computation of the generating 
series of the cumulants: 

TiX) = (AX + 1)(X + 1), 
^ ^ ^ ^ ^2(A + l 2n-l 

n>l 



[/(z) is the power expansion of i+[(a+i) +4A2j — _\ = q), where z i— > z-^/^ is the 

analytic version of the square root on the complement of the real non positive half line such that 

Theorem 3.5. Consider a G PiAp2 and define the formal power series M{X) = Y2n>i '>^2n{o,)X'^ 
and C{X) = "^^yi C2n(a)X". Then we have 

X 



1/2 



c = u 



{Xx{ToM)) 



<-i> 



Proof. By proposition 13.11 and the previous lemma, we have C o (Xx(To M)) = M. Hence 
roCo(Xx(ro M)) =ToM. DividingbyToM on both sides, one gets ^o(Xx{T o M)) = ^. 
Then by inversion, o (Xx(To M)) = X. At last, = (Xx(To M))'^^^'^ , which implies 

^ ° g = (Xx(ToM))<-^> ' then C = U ( (^,(^i,))<-i> - l) ■ □ 

3.4. Rectangular i?-transform of compactly supported symmetric measures. From 
now on, we will denote by (•)^^^ (resp. ^/■) the analytic version of the square root on the 
complementary of the real non positive (resp. nonnegative) half line such that 1^^^ = 1 (resp. 
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Consider a symmetric compactly supported probability measure /i. For n > 0, denote by 
m„(/x) its n-th moment and by C2n(/u) the number C2n{a) with a G PiAp2, where ^ is a (pi,p2)- 
probability space, such that for all A; > 1, ipi{{aa*)'^) = m2k{o)- 

Note that the generating series of the moments M^[X) = Yl,n>i '^2n{p)X'^ of p. is the power 
expansion of 

M^[z) = ^G^(^) - 1 {ze C\[0, +oo)), 

where is the Cauchy transform of fi on the upper half plane: G^{z) = J^^^^^^- Hence 
Xx(T o M) is the power expansion of 

2 

(z G C\[0,+cx))), 



H^{z) := zT o M^{z) = i^-^ + (1 - X)^zG^ (^-^^ 



that we shall call the rectangular Cauchy transform with ratio A of /x. Since ji is compactly 
supported, from the power expansion of the Cauchy transform G^, we know that is analytic 
in a neighborhood of zero, that Hfj_{0) = and that -ff^(O) = 1. It implies that we can invert 
the function H^^ in a neighborhood of zero and that the inverse function H^^ has a power series 
expansion given by (Xx [T o M^))^^^^. Note that i — 1 is analytic in a neighborhood of 

zero, with value zero in zero. At last, by theorem 13.51 Cfj_{X) := J2n>i '^2nia)X'^ is the power 
expansion of 



G^{z) := U ( — _^ — 1 ) (well defined for z small enough). 

Mil \z) 



where U{-) is the analytic function on a neighborhood of zero U{z) = +4Az] — 
if A = 0). 



Moreover, by (jlOp . we know that ^ G^ linearizes the rectangular convolution with ratio A: 
for all symmetric probability measures with compact support, 

(14) Gij,si^iy{z) = Gf,{z) + G„{z) {z smah enough). 

3.5. Rectangular i?-transform of unbounded measures: definition. The Cauchy trans- 
form and the rectangular Cauchy transform with ratio A of symmetric probability measures with 
unbounded support are defined with the same formulas as in the case of a compactly supported 
measures (see the previous paragraph). 

Proposition 3.6. Let A he a set of symmetric probability measures on the real line. Then the 
following assertions are equivalent 

(i) A is tight, 

(ii) for every < 9 < ir, lim ^H^{z) = 1 uniformly in fi & A, 

|arg2— ■7r|<9 

(iii) lim -Hnix) = 1 uniformly in fj, ^ A. 

x£{-oo,0) 

Proof, (i) ^ (ii) follows from the well known fact ( |BV93j ) that if A is tight, 
VO < ^ < TT, lim — G„( — ) = 1 uniformly in ri G A. 

z^O Z Z 
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{ii) =^ [iii) is clear. Suppose {in). Consider e > 0. Take rj > such that \/u € [0, 1], |An^ + (1 — 
X)u - 1\ < ri ^ \u - 1\ < e/2. Take x G {-oo,0) such that for all e A, |ii7^(x) -1\ <?]. Put 
-y/x = iy. Then for all £ A, 

l-G^(-) - 1| < e/2, i.e. / d/i(t) < e/2. 

Let M be such that Vt > M, > \. Then for all /x G A, /^([-M, M]^) < 2 -ftp^^^t) < e. 

□ ^ s/ 

Define, for a S (0, vr), /? > 0, A^^^ to be the set of complex numbers z such that | arg z — 7r| < a 
and |z| < /9. 

Let Ti be the set of functions / which are analytic in a domain Vj such that for all a G (0, vr), 
there exists [5 positive such that 

A family {fa)aeA of functions of W is said to be uniform if for all a £ (0, vr), there exists /3 
positive such that 

The following theorem has already been used, in other forms, to define Voiculescu's R- and 
^-transforms (see paragraph 5 of jBV93j). Its proof relies on Rouche theorem. 

Theorem 3.7. Let {Ha)aeA be a uniform family of functions ofTC such that for every a G (0,7r), 

lim — ^!^^ — I uniformly in a £ A. 
z^O z 
I arg 2— 7r|<a 



Then there exists a uniform family {Fa)a£A of functions ofH such that for every a G (0,7r) 

Fa(z) 

lim — = 1 uniformly in a £ A, 

z-*0 z 
I argz— 7r|<a 

and there exists f3 positive such that 

Va G A, HaO Fa = FaO Ha = Id on Aa,l3- 



Moreover, the family {Fa)aeA is unique in the following sense: if a family {Fa)aeA of functions 
ofTi satisfies the same conditions, then for all a G (0, vr), there exists (3 positive such that 



Va G A, Fa = Fa on A, 



For every symmetric probability measure on the real line ^, let us define the rectangular 
R-transform R^ with ratio A of /i by 

CJz) = u( 1 ) , 

where is defined by the previous theorem and the function U is the one defined at the end 
of section 13.41 
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One can summarize the different steps of the construction of the rectangular i2-transform 
with ratio A in the fohowing chain 



— 

sym. prob. Cauchy 
measure transf. 



CJz) = U 



H^{z) = XG^ (^) + (1 - X)V^G^ (^) 

auchy Vv^y Vv^y 



rect. /?-transf. with ratio A 



3.6. The special cases A = and A = 1. Note that the rectangular i?-transform with ratio 1 
(resp. 0), for a symmetric distribution /x, is linked to the Voiculescu transform of ^ by the 
relation C^{z) = y/zip ^{1 / y/z) (resp. C^(-z) = z(pp{l/ z), where p is the push-forward of /u by the 
function t ^ t^) (see paragraph 5 of |BV93] for the construction of the Voiculescu transform). 



3.7. Remark about the characterization of i?-transforms. The definition of the rectan- 
gular /^-transform of a probability measure ^ presents many analogies with the definition of its 
Voiculescu transform (/j^ (see [BV93j). So it seems natural to state, as the authors of [BV93j 
did for the Voiculescu transform (Proposition 5.6), a characterization of the functions that are 
the rectangular i?-transform of a symmetric probability measure. The characterization of the 
Voiculescu transform is based on the fact that for every probability measure the inverse of 
z + V^iiz) extends to a Pick function (i.e. an analytic function defined on the upper half plane, 
whose imaginary part does not take negative values), and on the fact that any Pick function 
equivalent to z at infinity is of the type 1/G^, so its inverse is z -|- ^^{z). But unless A = 
or 1, Pick functions do not appear in an analogous place in the definition of the rectangular 
i?-transform, so we cannot proceed similarly to characterize rectangular i?-transforms. 

Furthermore, the Levy Kinchine formula for a^^ (see jB-Glj ). compared with theorem 5.10 
of |BV93j ) will state that when is EB;^-infinitely divisible, there exists a unique symmetric 
ffl-infinitely divisible distribution v such that Cfj,{z) = ^/z(py(\l ^\ So the question of the 
characterization of rectangular iZ-transforms joins another question: can we extend the corre- 
spondence /i <-> to a bijective correspondence from the set of symmetric distributions into 
itself such that for all G^{z) = ^^fyiXj The analytic functions / on Aq,^^ of the type 
/(z) = \fzLpy(\l with V symmetric probability measure, are characterized by: 

(i) /(z) = 7(i), 

(ii) lim/(z) = 0, 

(iii) for all n and all z\^. . . ^Zn^ the matrix 

But nothing allows us to claim that the rectangular /^-transform of any symmetric distribution 
satisfies (iii), and that every function that satisfies (i), (ii), and (iii) is the rectangular R- 
transform of a symmetric distribution. 



/27(l + /(Zfc))-V^{l + /(2,)) 



is positive. 



kA = ^ 



3.8. Properties of the rectangular i?-transform. 

Theorem 3.8 (Injectivity of the rectangular i?-transform) . If the rectangular R-transforms with 
ratio A of two symmetric probability measures coincide on a neighborhood of in (— co,0), then 
the measures are equal. 
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Proof. If the rectangular i?-transforms with ratio A of two symmetric probabihty measures 
/i, coincide on a neighborhood of in (— oo,0), then by uniqueness of analytic continuation, 
they coincide on a A^^^, and so = Hi, on this set, and zH^(z) = zH^{z) on this set. So 
there exists M > such that for all y > M, 

X{iy2G^{iy)f + (1 - X)iyG^{iy) = X{iy2G,{iy)f + (1 - \)iyG,{iy). 

2 

But if p is a symmetric probability measure, iyG^{iy) = /f yi^t^ d/o(t) G (0,1]. So, by injec- 
tivity of It 1-^ An^ + (1 ~ -^)^ on (0, 1], G^{iy) = Gy{iy) for y > 1, and then, by analycity of the 
Cauchy transform, G^ = G^, and by injectivity of the Cauchy transform (see [XBlj or section 
3.1 of [HPOOj ) fj. = iy. □ 

The following remark gives a practical way to derive any symmetric probability measure on 
the real line fi from C^. 

Remark 3.9 (How to compute fi when we know ?). First, we have z/H^^{z) = T{Ci^i{z)), 
for z G C\M^ small enough. From this, we can compute H^{z) for z £ C\M^ small enough. 
Then we can use the equation, for z £ C\M"*', 

li^,W = A(-LG,(-L))%(i-A)-Lo„(-L). 



Moreover, when z £ C\M''~ is small enough, l/\^z is large and in C , so -^G^ 
to 1. ^Hf^{z) is also closed to 1, and for h,g complex numbers closed to 1, 



is closed 



h = Xg^ + {l-X)g^g = V{h), with V{z) = ^ ^ +4Az)2 ^ _ i) + i. 

So one has, for z G C\M^ .small enough, -^G^{-^) = ^(— which allows to recover /i 

The following lemma is an easy consequence of proposition 13.61 and theorem 13.71 

Lemma 3.10 (Tightness and rectangular i?-transform) . Let A be a set of symmetric probability 
measures. Then we have equivalence between 

(i) A is tight, 

(ii) for any < a < vr, lim C^{z) = uniformly in fj, £ A, 

|arg z—iT\<a 

(iii) lim Ra{x) = uniformly in fi £ A. 
x—*0 

xe{-oo,0) 

Theorem 3.11 (Paul Levy's theorem for rectangular i?-transform) . Let (^„) be a sequence of 
symmetric probability measures. Then we have equivalence between: 

(i) (/x„) converges weakly to a symmetric probability measure; 

(ii) there exists a, (3 such that 

(a) lim Ca^{z) = uniformly in n, 
z—*0 

I arg z—Tr\ <a 

(b) the sequence (C^„) converges uniformly on every compact set o/A^^^^; 

(iii) (a) lim Cu^(x) = uniformly in n, 

X — ^0 
a::e(-oo,0) 

(b) there exists (3 > such that the sequence {C^^) converges pointwise on [—(3,0). 
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Moreover, in this case, denoting by /x the weak limit of (nn), for every a, there exists (3 such 
that the sequence (C^,J converges uniformly to on every compact set of Aa^p- 

Proof. (z) =^ (ii): suppose that converges weakly to /i. Then by the previous lemma, 
we have (a) of (ii). So there exists (3 > such that on Aa,/3, for all n, |C^„| < 1. So, by 
Montel's theorem, it suffices to show that the only possible limit, for uniform convergence on 
every compact, of any subsequence of (C^„) is C^. Let C be an analytic function on such 
that a subsequence (C^^ ) of (C^„) converges uniformly to C on every compact of Aq,^^. We 
know (see |A61j or section 3.1 of [HPOOJ) that the sequence (G^^) converges uniformly on every 
compact of the upper half plane to G^. So the sequence converges uniformly on every 

compact of the complement of [0, +cxd) to H^- Since limC(2:) = 0, to prove C = Cn, it suffices 

to prove that 



(AC+1)(C + 1) 
But (AC^(z) + 1) + 1) 



(AC^ + 1) (C^ + 1) . 



So it suffices to prove that 



{XC{z) + 1) + 1) 
so by theorem [331 it suffices to prove that 

z 



We have 



(AC(z) + 1) {C{z) + 1) 



z. 



(AC(^)+l)(C(^)+l) 



(AC(^)+1)(C(^)+1) 



< 



+ 



{XC{z) + 1) (C(z) + 1) 



{XC,,Jz) + l) {C,,Jz) + l) 



(1) 



{XC,,Jz) + l) {C,,Jz) + l) 



H 



{XC,,Jz) + l) {C,,Jz) + l) 



(2) 

By continuity of H^, (1) tends to zero when n tends to infinity, and, since Hf^^ converges 



uniformly to i?^ on every compact, (2) tends to zero. So if^ \^ (AC(z)+i)(C(^)+i) 

(ii) ^ (iii) is clear. 

(iii) ^ (i): suppose (Hi). Then by (i) (ii), every limit of a subsequence of (/x„) has a 
rectangular i?-transform equal to the pointwise limit of (C^„) on [/?, 0). By uniqueness of analytic 
continuation, all the limits of subsequences of (/U„) have the same rectangular i?-transform, so, 
by injectivity of this transform (previous theorem), there cannot be more than one limit of 
subsequence of (/^n)- As by (a) of (iii), the set {fin ; n G N} is tight, (/x„) converges weakly to 
a symmetric probability measure. □ 
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3.9. Rectangular convolution of measures with unbounded support. 

Theorem 3.12. The binary operation defined on the set of compactly supported symmetric 
probability measures in section\^ extends in a unique way to a commutative, associative, and 
continuous (with respect to the weak convergence) binary operation on the set of symmetric 
probability measures on the real line. This operation is called the rectangular free convolution 
with ratio A. Moreover, for all n, v symmetric probability measures one has C^^^u = Cf^ + C^,. 

Proof. Let /U, u be symmetric probability measures. If (resp. (fn)) is a sequence of 
compactly supported symmetric probability measures that converges weakly to /i (resp. z^), 
then by theorem 13.111 and ([HI), the sequence (nn^^^n) converges to a measure whose rectangular 
i?-transform is + (thus by injectivity of the rectangular i2-transform, this measure does 
not depend on the choice of the sequences (/i„) and (vn) and is equal to na^u when fi and ly 
are compactly supported). Note that C^^^u = + C^, stays true for all symmetric probability 
measures ^ and i^. Moreover, this equation shows that is a commutative, associative, and 
continuous binary operation on the set of symmetric probability measures. □ 

Now we can extend the proposition 12 . II to the case where fii, H2 are not compactly supported: 

Theorem 3.13. Consider sequences qi{n),q2{n) such that 

qi{n) <q2{n), qi{n) — > +oo, qi{n)/q2{n) — > Ae[0, 1]. 

Let, for alln, An,Bn be independent qi{n)xq2{n) random matrices, one of them being bi-unitarily 
invariant, such that the symmetrizations of the singular laws of An,Bn converge in probability 
to probability measures /i, Then the symmetrization of the singular law of An + Bn converges 
in probability to nm^v. 

To prove the theorem, we will need the following lemma, for which we need to introduce 
an extended functional calculus. When F is a real Borel function on the real line, for any 
qi{n)xq2{n) matrix C, we define F{C) to be the matrix 

l<i<iJ2(n) 

where /ii, . . . , > 0, and f7, F are respectively qi{n)xqi{n), q2{n)xq2{n) unitary matrices 
such that 

l<J<g2(n) 

Note that {F + G){C) = F{C) + G{C) and F{C)G{C)* = when FG = 0. For any probability 
measure a, we will denote by F{a) (for example a'fi, \a\.. . ) the push- forward of o" by F. 
We denote by fin the spectral measure of an hermitian matrix H. Recall that we denote by 
a the symmetrization of any probability measure a. For example, the symmetrization of the 
singular law of a rectangular matrix A will be denoted by I1\a\- 

Lemma 3.14. Let (Gn) be a sequence of qi{n)xq2{n) random matrices such that jl\c,^\ converges 
in probability to a probability measure cr. Then, for any odd real function F continuous at 
a-almost every point o/M and such that F{M~^) C M"*", fl^p(^(j^'^^ converges in probability to F{a). 

Proof. F{a) is symmetric and the symmetrization is a continuous operation, so it suffices 
to prove that fJ-\F{c„)\ converges in probability to \F{a)\ = F{\a\). But M|f(c„)| = -^(/^|C„|)> 
fJ-\Cn\ converges in probability to \a\, so the conclusion follows from the fact that the function 
a I— > F{a) on the set of probability measures on the real line is weakly continuous at \a\ because 
F is continuous at |(7|-almost every point of the real line ( [B68] ). □ 
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Proof of the theorem. First of all, for Un, Vn independent gi(n)x q2{n)xq2{n) uni- 
form random unitary matrices which are independent from the random matrices 
UnAnVn,UnBnVn are independent (it is a property of the Haar measure on compact groups), 
have the same singular values as respectively An,Bn, both of them are bi-unitarily invariant 
and the singular values of UnAnVn + UnBnVn = Un{An + Bn)Vn are the same ones as the ones 
of An + Bn. Hence one can suppose that both An, B^ are bi-unitarily invariant. 

The function that maps a probability measure a on to the symmetrization a^/^ of its 
push-forward by x — > x^^ is continuous. So it suffices to prove that the push-forward of the 
spectral law of \An + Bn\ by x — > x^, i.e. the spectral law of M„ := {An + Bn){An + Bn)* 
converges in probability to (fim^^v)'^ . We can define a distance on the set of probability measures 
on the real line with the Cauchy transform by (cri,(T2) sup{\Gai{z) — Ga2{z)\ ; > 1}. 
This distance defines the topology of weak convergence ( |A61] . |PL02j ). The Cauchy transform 
of the spectral distribution of an hermitian matrix M is the normalized trace of its resolvant 
9^^(M) = (z- M)~^ So the spectral distribution of a sequence (Xn) of qi{n)xqi{n) hermitian 
random matrices converges in probability to a probability measure a on the real line if and only 
if for each e > 0, the probability of the event 

sup |tr9^^(X„) -G^(z)| > e 

Sz>l 

tends to zero as n tends to infinity. Choose e > 0. We will show that 
(15) P{sup|tr(^n,(M„))-G(^a,.)2(^)| >e}™0, 

where P designs the probability measure of the probability space where the random matrices 
are defined. Let us define, for t > 0, the function Ft on the real line by 



Ft : X 



X if |x| < t, 
else. 



For every probability measure cr, Ft{(j) converges weakly to a when t tends to infinity, so, by 
continuity of m^, there exists t € (0, -|-oo) such that t and —t are not atoms of the measures /i 
and and such that 



(16) sup |G(i.^(^)a^j.,(^))2(z) - G(^H^^)2(z)| < 

(17) ;u(M\[-t-l,t + l]) + i/(]R\[-t-l,t + l]) < 



e 
3' 

e 

18' 



We win now use the notations Ft = F, G{x) = x - F{x) and Mn,t = iF{An) + F{Bn))iF{An) + 
F{Bn))* ■ By triangular inequality and (fT6l) . we have 

P{sup|tr(in,(M„))-G(^B,.)2(^)| > e} 

S>2>1 

< P{sup|tr[9^,(M„)-9^,(M„,t)]| > §} 

S2>1 

(18) + P{sup|tr9^,(M„,t) -G(^(^)ffl^^(,))2(z)| > §} 

To treat the first term of the right hand side of (|18p . recall that for any qi[n)xqi{n) matrix T, 
|trT| < 1 1^1 1 rg T. The operator norm ||5H2(M„) — y{z{^n,t)\\ is not less than 2, because 
9z > 1. Moreover, due to the equation 

1H,(M„) - lH,(M„,t) = -<K,(M„)(M„ - M„,t)in,(M„,t), 
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the rank of yiz{Mn) — ^ziMn^t) is not more than the one of M„ — Mn^t- One has (omitting the 



indices n in An,Bn) 

Mn - Mn,t = G{A)G{A)* + G{B)G{B)* + G{A)B* + F{A)G{B)* + G{B)A* + F{B)G{A)*, 
so its rank is not more that 3vgG{An) + 3rgG(i3„). Hence 

P{SUP \tl[9\z{Mn)-9\z{Mn,t)]\ > §} 

S2>1 

< p{-|-^(rgG(A„)+rgG(S„))>|} 

= F {f,AAm-t,t]) + f^BAm-t,t]) > 

which tends to zero when n goes to infinity, by (I17p . On the other side, the second term of right 
hand side of ()18p goes to zero when n goes to infinity by definition of and by the previous 
lemma. So (llSp is proved. □ 

3.10. Examples. In this section, we give examples of rectangular /^-transforms of symmetric 
probability measures, and examples of computations of rectangular free convolutions. 

The case where A = 1 (see section [321), where the rectangular i?-transform and the rectangular 
free convolution are essentially the Voiculescu transform and the free convolution defined by 
Voiculescu and Bercovici, doesn't present any new interest (one can find the little number of 
examples known in, e.g., [HPOOj ) . so we will suppose that A € [0, 1). Unfortunately, for A G (0, 1), 
few computations can be done, still less than for the well known case A = 1. 

3.10.1. Convolution of symmetric Bernouilli distributions. 
Proposition 3.15. Suppose A > 0. Then r := support 
(19) [-(2 + -(2 - u [(2 - k)^/\ (2 + k)^/% 

with K = 2(A(2 — A))-*^/^ G (0, 2), and it admits a density with formula 

[^2 _ (^2 _ 2)2] 1/2 

7rA|x|(4-x2) 

on its support. 

Remark 3.16. By theorem \3.13\ this proposition means concretely that for sequences qi{n) < 
q2{n) which tend to infinity such that qi{n)/q2{n) tends to X G (0, 1), if one considers independent 
random matrices An, Bn,Un,Vn such that for all n, An,Bn are qi{n)xq2{n) matrices such that 
AnA^, BnB^ are the identity matrix and Un, Vn (^re uniformly distributed unitary random 
matrices with respective sizes qi{n),q2{n), then the symmetrized singular law of An + UnBnVn 
tends to the probability measure r with the density of (f20l) on the support of (fT9]l . Equivalently, 
one can say that the spectral law of 

(An + UnBnVn){An + C/„5„F„)* 

tends to the push forward of t by the map x i— > x'^ , i.e. to the probability measure with support 
[2 — K, 2 + k] and with density 

[^2_(^_2)2]l/2 

7rA2;(4 — x) 



RECTANGULAR RANDOM MATRICES, RELATED CONVOLUTION 25 

Remark 3.17. Note that the rectangular free convolution of the Dirac mass at zero with any 
other symmetric probability measure is the measure itself, and that for any c > 0, '^-e+'?e ^^ <^-c+'^e 

can be deduced from this proposition using a dilation. Note also that if X = 0, then r = ^-^^^^ . 
Proof. Define i/ = We have G^{z) = hence G^{1/^/^) = ^ and 

^-^') - (13^ + ^1 - ^)JT^ - (l-z)2 • 

Hence for x,y £ C\R^ in a neighborhood of zero, the equation H^{y) = x is equivalent to 
(A — + y = X — 2xy + xy^, i.e. to 7/^ (x + 1 — A) — y(l + 2x) + x = 0. In this case, since 
y = ifj7^(x) must tend to zero as x tends to zero, one has 

1 + 2x - (1 + 4Ax)i , X , (1 + 4Ax)^ - 1 + 2x 

H,. (x) = ^ r and 1 = . 

^ ^ 2(x + l-A) Hu\x) 2 



Note that (A + l)^ + 4A - 1 j = (^A + (1 + 4Ax) 2 j , hence 



-A-1 + 
C,{x) = 



(A + 1)2 + 4A 



1 



2A 2A L 



:i + 4Ax)3 - 1 



Thus by (jl4p and remark [3T9l for y G C\M''" in the neighborhood of zero, if one denotes (l+4Ay) 
by w, one has 

Xy {w^ - l)/4 



l + 4Ay + (A-l)(l + 4A2/)3 u;2 + (A-l)u'' 



1 



Hence for x G C\M+ in the neighborhood of zero, if one denotes Ht-{x) by y and (1 + 4Ay)2 by 
w, one has H~^{y) = x, i.e. w"^ — 1 = Axw^ + 4x(A — i.e., since w tends to 1 as x tends to 
zero, 

w = ^^^[2x(A - 1) + [1 - 4x + 4(A - 1)^x2]!]. 
Then an easy computation leads to 

(A -1)2 + 1^= / A-l + 2[l-4x + 4(A-l)V]i y 

Thus by remark [3.91 for z G C\M^ in the neighborhood of zero, 

1 (—\ = V (^^\ ^ ^ A _ ^ ^ A - 1 + 2[1 - 4z + 4(A - ifz^] h \ 



z \^/z \ z 2X\ 1 - 4z 



, 1 



hence ( ^) = -V A - 1 + A-i+2fi-4/(i/yy +4(A-i)V(i/^)^]^ , ^ ,1^- ^^^g^ 

^ ^ . ^ (2^^-4)(A-l)+2[(^2-2)2-,.2]^^' f ^ ^- 

enough, GT-(^j = 2\z{z^-A) — ^ ^ analytic contmuation, for any z G (L , 

Gr{z) = ^^~'^\ll^^^}lt^'^^^'^ 1 where 5(2:) is analytic on C~ and satisfies 5(2:)^ = {z"^ — 2)^ — k? . 
In order to compute r, we shall use the following well known result (lemma 2.17 of |B06] ): for 
Lebesgue-almost all real number x, the limit, as z = n + if G C~ tends to x in such a way that 
{u — x)/v stays bounded, of ^Gr{z)/T: is the density of the absolutely continuous part of r at x. 
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Let us define 

O := (-(2 + k)^I^, -(2 - U ((2 - k)^/^, (2 + k)^'^). 

By what precedes, for Lebesgue- almost all x G M, there is no sequence Wn > tending to 
such that lim QGr{x — iwn) < 0. Hence since for all x £ O, (x"^ — 2)-^ — is negative, for 

n— >oo 

Lebesgue-almost all x G O, there is e > such that for all w £ {0,e), for z = x — iw, 



( (4-2^^)(A-1)+2V(z2-2)2-k2 ^„ 
2A2(4-^2) n J, > U, 



(4-222)(A-1)-2V(22-2)2-k2 -r^r. 

2A2(4-22) n X ^ u. 



It implies that for Lebesgue-almost all x G O, 

c^^.^.w 5(V(x^-2)^-K^) [k^-(x^-2)^]^/^ 

Now, to complete the proof, it suffices to prove that r is equal the absolutely continuous part 
of its restriction to O, i.e. that 



o 



[^2 _ (^2 _ 2)21 1/2 

7rA|x|(4 - 



-dx = 1. 



This can be proved using the changing of variable y = -r^ — o\2\i/2 the well known formula 

(K —(X —I) ) 

of the Cauchy transform of the Cauchy distribution. □ 

3.10.2. The rectangular analogues of the Gaussian distributions. Here, A > (the case A = 
can be found in |B-Glj ). It is proved in [B-Gl] that for all cr > 0, the rectangular analogue of 
the symmetric Gaussian distribution with variance o"^ is the law N^2 with density 

[4A-(a^x^-l-A)^]^/^ , 

2^^A^^N ""^^ " ^' 

where x stands for the characteristic function of the interval [{1 — X^^'^)'^ , (1 + A^/^)^]. Its support 
is [-cr(l + X^/^),-a{l - A^/2)] y [^(x _ Ai/2),o-(l + X^^"^)], its rectangular i?-transform is a^z. 
One has, for all a,a > 0, N^2m^N^2 = N^2j^^2. 

Remark 3.18 (Growth of the support in the related semigroup). Note that in the additive 
semigroup {{Nc ; c > 0}, m^) the size of the support of Nc in not linear in c but in C2 . This kind 
of phenomenon, which recalls us that we are in a gaussian context, but still surprising since the 
support (and not the variance) is concerned, had already been observed in free probability. 

3.10.3. The rectangular analogues of the Cauchy distributions. For all t > 0, the rectangular 
analogue of the symmetric Cauchy law with parameter t by the Bercovici-Pata bijection with 
ratio A is 



t 

^* t(i~A) ,(x)- 



\_t'{x-ir' 



4x2 



2 

dx. 



Its support is M\(— ^^^2"^! ^^^y^). It is proved in |B-Glj that the i2-transform with ratio A of 
ft admits an analytic extension to C\M^ given by the formula it^fz. Hence one has, for all 
t\M > 0, iyt-^m^iyt2 = fti+ta- 
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3.10.4. The rectangular analogues of other symmetric stable laws. It can also be proved that for 
all a € (0, 2), the rectangular analogue of the classical semigroup of stable laws with index a is 
the set of laws with rectangular i?-transforms with ratio A given by the formula z — C(— z)"/^, 
with C > and where the a/2-th power is defined with respect to the canonical determination 
of the argument on C\M~. Unfortunately, though we give their rectangular ii-transforms with 
ratio A, except for a = 1 (see the previous paragraph), we cannot give the densities of these 
laws. However, just like for free stable laws (see [BPB99] ). results can be proved about their 
densities and supports (see |BBG07j ). 

3.10.5. The rectangular analogues of the symmetric Poisson distributions. It is proved, in [B-Glj . 
that for all c > 0, the sequence 

has a weak limit, that we shall denote by Pc, and that the rectangular i?-transform of Pc is 
given by the formula We cannot compute Pc itself, but we know that for all c, c' > 0, 

Pc^\Pc' — Pc+c' ■ 



4. Appendix: a technical result 

In this section, we do not make any particular hypothesis about pk^s: they are nonnegative, 
and at most one of them is zero. So we can suppose that p2, ■ ■ ■ , Pd ai^e positive. 

For all integers n,k, for each partition V of [k], [n]^ will design the set of elements i = 

(ii, . . . ,ik) of [n]'^ such that h Z h' <^ i^ = ih'. Moreover, \V\ denotes the number of blocks of 
V. 

Definition 4.1. An nxn matrix is said to be an almost diagonal matrix if it is of the type M , 
with M diagonal (either square or rectangular) matrix. In this case, for all j £ [n], Aj will 
denote the only nonzero term on the j-th column of A if it exists, and zero in the other case. 

Lemma 4.2. Let K be a positive integer, and let for each n > 1, A{l,n), . . . , A(K,n) be a 
family of almost diagonal nxn matrices. Let TC be a subset of [K], and let V be a partition of 
[K] such that for all {h} singleton class ofV, A{h,n) is a diagonal matrix and TrA{h,n) = 0. 
We suppose too that the family {A{i,n))i^n is uniformly bounded, and that for all h €TC, for all 
i G {qiin) + 1, . . . , n}, Ai{h, n) = 0. Then 

K _ 
llAi^{h,n) = 0(9i(n)l^«l-P«/2„l^«l-P«/2)^ 

where Vn = {B nTi; B e V}\{f/>} an_d Vn{B n {K\n) ; B G V}\{f/>}) and pn (resp. p^j the 
number of singletons in Vn (resp. ofVn). 

Proof. We will prove this result by induction on p = p-^i + p^. If p = 0, the result is clear. 
Now, let us suppose the result to be true to the ranks 0, ... ,p — 1, with p>l. Consider V such 
that Pn +P'H — P-> ™d a singleton class {/iq} of V. Denote by V' the partition V — {{h^}} of 
{1, . . . , /iQ, . . . , and let, for each class B olV with B ^ {/io}i ^(^o ~^ B) be the partition of 
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[K] obtained from V by linking the classes {ho} and B. We have 

( ^ ^ 

J2 Ai{ho,n) Yl Ai^{h,n) 

Yie[n]-{ii,...,i^Q,...,ix} J 



h=l 



But, since A{ho,n) is diagonal and has null trace, 

Y Ai{ho,n) = - Y Ai{ho,n). 

ieln]-{ii,...,ihg,--;iK} ie{ii,—,iho,--;iK} 
So ^ ^ 

l[A,,{h,n) = -Y E Il^Hih,n). 

For example, suppose ho £ 7i (the other case is treated in the same way). With the induction 
hypothesis, and dividing the sum in i? G Vn and B G Vn^ one has (each es being 1 or 
according to whether B is a singleton or not) 

Y llAi,{h,n) = Y 0(gi(n)l^«l-i-(^'«-^-^«)/2nl^«|-J^«/2) 

ie[n]'P/i=i BeVn 

B^{ho} 

+ Y O (^gi(n)l^«l-(P«-^)/2nl^«|-^-(^«-"-B)/2) 
Bev-H 

= qi{n)-^ Y O (gi(n)l''«l-^'«/2nl^«l-^«/2) 

BeVn 
B^{ho} 

+gi(n)^n^ Y O (^qi(n)^^»^-P^V»^-Pnf^^ . 
BeVn 

Since qi{n) = 0{n), the result is proved. □ 

The following proposition, main result of this subsection, is divided in two cases, which give 
very similar results with very similar proofs. The proposition involves special classes of matrices. 

Definition 4.3. For k G [d], we define Ufe(n) to be the set ofnxn matrices whose {k,k)-th block 
is a qk{n)xqk{n) unitary matrix, and with other blocks zero. 

Note that (l[Jfc(n),x) is a compact group, isomorphic to the one of qk{n)x qk{n) unitary 
matrices, and random matrices of Ufc(ra) distributed according to the Haar measure will be said 
to be uniform. In the following, we are going to use an integer A'^ and the cyclic order on [N]. 
This means that to put the index A?^ + 1 on an element is equivalent to put the index 1. 

Proposition 4.4. Let, forn > 1, V{s,k,n) (s GN,k G [d]), be a family of independent random 
matrices, such that for all s, k, V{s, k, n) is uniform on Ufc(n). Fix N positive integer and R> 0. 

Case 1°) Let, for each n E N, D{l,n), . . . ,D{N,n) be nxn constant matrices such that for 
all r € [A^], ||D(r,n)|| < R, and there exists u{r),v{r) G [d] such that one of the two following 

conditions is realized: 

(i) u{r) = v{r) and for all n, D(r, n) = Pu(r)('^)j 
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(ii) for all n, D{r,n) is an almost diagonal matrix of the type M, where M is a qu{r){''^)^ 
qy{r){n) diagonal matrix, andE{D{r,n)) — > 0. 

Consider (si, A;;);g[7v] G (Nx[(i])^, and suppose that for all r £ [N] such that (i) above is satisfied, 
{sr, kr) 7^ (sr+i) ^r+i)- Consider also mi, . . . ,mj\f G Z — {0} and r] > 0. Then the probability of 
the event 

{\\E{V{si,ki,nr'Dil,n) ■ ■ ■ DiN - l,n)VisN,kN,nr'' D{N,n))\\ <r,} 

tends to 1 as n goes to infinity. 

Case 2°) Let, for each n £ N, D{0, n), . . . , D{N, n) be nxn constant matrices with norms also 
< R. We suppose that there exists v{0),u{N) G [d] such that D{0,n) (resp. D{N,n)) is an 
almost diagonal matrix of the type M, where M is a qi{n)xq^((^-^{n) (resp. '?u(Ar)(ra)xqi(n)j, and 
that for all r £ [A'" — 1], there exists u{r),v{r) £ [d] such that one of the two following conditions 
is realized: 

(i) n(r) = v{r) and for all n, D{r,n) =Pu{r){''^)> 

(ii) for all n, D{r,n) is an almost diagonal matrix of the type M, where M is a qu{r)i''^)^ 
qv(r)in') diagonal matrix, andE{D{r,n)) — > 0. 

Consider (sj, fci)ze[Af] ^ ['^D^' '^^^ suppose that for all r £ [A^ — 1] such that (i) above is 
satisfied, {sr,kr) ^ {sr+i,kr+i). Consider also mi, . . . ,m]\f G Z — {0} and r] > 0. Then the 
probability of the event 

{II E(L>(0, n)V{si, ki,nr'D{l, n) ■ ■ ■ D^N - 1, n)y(sjv, kN,nr^D{N, n))|| < rj} 

tends to 1 as n goes to infinity. 

Remark 4.5. By linearity ofE, the result stays true in case 1°), (resp. in case 2°)) if for all 
r in [N] (resp. in [N — 1]) such that (ii) is satisfied, B{r,n) is replaced by a sum of constant 
almost diagonal matrices whose images by E also tend to zero. 

The proof is inspired from the one of Theorem 4.3.1 p. 147 in [HPOOj . 
Proof. Both proofs wiU be made together, we wih only have to separate them sometimes. 

Step I. Cases 1°) and 2°) Let us denote, for all r which satisfies (ii) above and u{r) = v{r), 
D{r,n) = D'(r,n) + Xr{n)pu(r){n'), with E(L>'(r, n)) = 0. 
Then by hypothesis, \h{n) — > 0. Hence, by linearity and boundedness hypothesis, it suffices 

n— >oo 

to prove the result for D{r,n) replaced by D'{r,n). Thus we can from now on suppose that for 
all r such that (ii) is realized, we have E{D{r,n)) = for all n. 

Step II. Case 1°) By definition of E, it suffices to prove that for all k £ [d], the normalized 
trace of the k-th diagonal block of the product 

V{si, ki,n)'^'D{l, n) • • • D{N - 1, n)V{sN, kN,n)'^^D{N, n) 

converges in probability to zero, i.e. that 

(21) ^ Tipk{n)V{si, ki,nr^D{l, n) ■ ■ ■ D{N - 1, n)V{sN, kN,nr^D{N, n)pk{n) 
converges in probability to zero. So let us fix /c G [d], and let us prove it. 
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li ki ^ k or v{N) ^ k, then the A:-th block of the matrix is always zero, so the result is clear. 
So let us suppose that ki = v{N) = k. We can then remove pfc(n) in (|2T|) . By Markov inequality, 
it suffices to prove that ()2ip tends to zero in L^. Since for all matrix M, | Tr = Tr TV M* , 
we only have to prove that the expectation of 
(22) 

Tr V{si,ki,n)"'^D{l, n) ■ ■ ■ D{N, n) Tr D{N, n)*V{sN, kN,n)-"'^D{N - 1, n)* • • • V{si, ki,n)-"'^ 

is o{qk{nf). 

Case 2°) We have to prove that 

Tr D(0, n)V(si, ki,n)"''D(l, n) ■ ■ ■ D(N - 1, n)V(sN, kN,n)'^^D(N, n) 

converges in probability to zero. Since Tt XY = TiYX, it suffices to prove that 

(23) ^— Tr V(si,ki,n)'^^D(l, n) ■ ■ ■ D(N - l,n)V(sN, kN, n)'^^D(N, n)D(0, n) 

converges in probability to zero. From now on, we will denote D{N, n)D{0, n) by D{N, n), which 
is now an nxn matrix of the type M with M qu{N)i''T')^Qv{o)i''T')j but which has not more than 
gi(n) nonzero entries. For the same reason as above, we only have to prove that the expectation 
of 
(24) 

Tr V{si,ki,n)"'^D{l, n) ■ ■ ■ D{N, n) Tr D{N, n)*V{sN, kN,n)~'^^ D{N-l, n)* ■ ■ ■ F(si, /fci, n)~'^i 
is o(gi(n)2). 

Step III. Cases 1°) and 2°) Let us expand both traces in the previous product, and re index 
the sum on partitions. 

Let M = \mi\ H h \mN\, let M(0) = and, for r £ {1,. . . ,N}, 

M{r) = \mi\ + • • • + \mr\, 
let M{N + 1) = 2M, and let, for r G {iV + 2, . . . , 2N}, 

M{r) = M + |mAr| + \mN-l\ H h \m2N+2~r\- 

Let also, for r G {iV + 1, . . . , 2N}, D{r, n) = D{2N + l-r,n)*. We write, for k G Nx [d], 

Vij{K,l,n) = V{K,n)ij, 
Vij{K,-l,n) = V{K,n)j^i. 

With those notations, there exists two functions k and e such that the expectation of (I22p is 



2N 



ie[n]2« r=l 



2M 



.h=l 



where E denotes expectation and for all i = (ii, . . . , i2M) G [n 



\2M 



jl = j2 = h, ■ ■ ■ ijM = h, JM+I = iM+2, jM+2 = «M+3, • • • ,j2M-l = «2M, J2A/ = ^M+l- 
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Let us introduce the partition Q of [2M] defined by /i ~ /i' if and only if if^ = i^'- Then we 
can rewrite the preceding sum 



2N 



Q partition ig[nl2 r=l 
of [2M] 



2M 



Wvi^,3h{i^{h),e{h),n) 



.h=l 



Thus it suffices to prove that for ah partition Q of [2M], the sum 



2N 



2M 



(25) E n^(^'^)^MM^ 

ie[n]S r=l L?t=l 

is o{qk{n)'^) in case 1°), and o{qi{rif) in case 2°), as n goes to infinity. 
So we fix a partition Q of [2M]. 

S'tep IV. Cases 1°) and 2°) Let us denote, for all h G [2M], = {s{h),k{h)). Note that, 
for i G for 

(26) E 



2M 



n ^'hjh('^(^)'^(^)''^) 



./i=i 

to be nonzero, i has to satisfy 

(27) V/i G [2M],ih,jh e {9i(?^) H h qk{h)-i{n) + 1, • • • ,gi(n) H h qk{h)~i{n) + qk(h){n)]. 

Indeed, in the other case, one of the factors in the product is zero, by definition of M i— > M. 
Note also that for all elements i of [n]®, which satisfy (|27p. the expectation (j26|) is the same, by 
invariance of Haar measure on unitary groups under permutation of rows and columns. Let us 
denote this common expectation by Eg. So, one can write the sum of (j25p 

(28) ¥.qY,X{D{ 

i r=l 

where the sum is taken on elements i of [n]® which satisfy (j27p . 

It is well known (see |HPOO| . equation (4.2.11)) that if, for each m, U{m) is a uniform unitary 
771 X ui random matrix, and 1 ^ jm < fn, then the sequence ^{\U {m)i^jj\'^^') does not 
depend on the choices of im,jm and 

(29) E(|C/(m),„,,„|2*^)^^*^ = 0(m-i^). 
So, by Holder inequality, for i G [n]^, 

2M 

IEqI < llnK,jMh),e{h),n)r)'^'-' = 0{q,{n)-^'' n-^'"), 
h=l 

where M' = |mr| and M" = \mr\. Thus it suffices to prove that 



l<r<Ar 

kr = l 



l<r<N 



2N 



(30) 



n 



jM(r) ' 



i r=l 



is o{qk{n) qi{n) n ) in case 1°), and o{qi{n) ~^ n ) in case 2°), where the sum is taken 
on the elements i of [n]^ which satisfy ([27]). 
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Step V. Cases 1°) and 2°) We are going to prove it as an application of lemma 1121 The sum 



of (j30j) is the sum, on the same i's, of 

D{1, n),,,,^,^^ D{2, n),^,,,,^^ ■ ■ ■ D{N - I, n),,,,^_,,^,Z)(iV, n)^, 
X DiN + 1, n)i,,^,Z)(Af + 2, n)i^^^^^^^DiN + 3, n)i^^^^^^ ■ ■ ■ D{2N, n)^^,^,^, . 



i.e. of 

*M(1) + 1 -^(^' '^)«M(2) + l ' ' ' D(N 1) '^jjjvffiv-lj + l 

X DiN + l,n)i,,^,DiN + 2, n),,,(^^,,I)(iV + 3, n),,,,^^3, • • • Z?(2iV, n)^,,^,^, . 
Adding some matrices of the type Pk{n) in the product, one sees that this sum is equal to 

2M 

(31) E n^*'^(^'")' 

ig[n]Q h=l 

where for /i e [2M], 

' Pk{h~i){n)D{r,n)pk(^h){n) if 3r G {0, . . . , - 1}, 

= M(r) + 1 

Pfe(h)H if /i G [M] and /i ^ 

{M(0) + 1,...,M(A-1) + 1} 
Pfe(2M) (n)L'(A + 1, n)pfc(M+i) (n) if /i = M + 1 

if 3r G {A + 2,...,2A}, 

h = M{r) 
if /i G [2A] and h ^ 

[M + 1] U {M{N + 2), . . . , M(2A)} 



A{h,n) 



Pk{h)in)D{r, n)pfc(^+i)(n) 



Now, note that by lemma 4.2.2 of jHPOOj (or because the uniform distribution on the unitary 
group is left and right invariant), if i G [n]^ corresponds to a nonzero term in equation (I25p . 
then for all k G Nx [d], for all 1 < q, /? < n, 

\{h; ih = a,K(h) = K,e(h) = l}\ = \{h; jh = a,K{h) = K,e{h) = -1}\, 
\{h; jh = P,K{h) = K,e{h) = l}\ = \{h; ih = P,K{h) = K,e{h) = 

Thus, if {h} is a singleton class in Q, and if Eg is non null, one has, using the cyclic orders on 
[M] and on {M + 1, M + 2, . . . , 2Af}, 

e{h) = 1 K{h) = K{h - 1), e(/i - 1) = -1, 
and e{h) = -1 k(/i) = - 1), e(/i - 1) = 1. 

It clearly follows that there exists r G {0, . . . , 2A} — {A^ + 1} such that h = M{r) + 1. Moreover, 
by hypothesis, in case 1°), A{h,n) = D{r,n) is a diagonal matrix with null trace, and, in case 
2°), A{h, n) = D{r, n) is a diagonal matrix with null trace whenever h ^ {M, M + 1}. 

Case 1°) Let us apply the lemma with V = Q and TC = {h £ [2M] ; k{h) = 1}. By definition 
of the A{h, n)'s, for all h G 7i, for all i G {qi{n) + 1, . . . , n}, Ai{h, n) = 0. By the lemma, the sum 
of ([30]), which is equal to the sum of is 0(gi(n)l^«l-P«/2j^|PH|-PH/2). But for any partition 
X with X singletons of a set 5, one has \X\ < x + {\S\ — x)/2, so 1^"! — x/2 < \S\/2. So, 

\'Pn\ + \rn\-^-^< {2M)/2 = M' + M" , 
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Moreover, since for all h G [2M], k{h) = 1 implies /i G \V-h\-PhI'^ - ^ ''^^ ' ^^^"^ ^ ^il = ^" ■ 
Thus, since qi{n) < n, nl^^l-P^/^ = „FHhP«/2-M"^M" ^ 0((7i(n)l^«l-PH/2-A^"n^^"). Hence 

The sum of pop is 0(gi(n)*^'n*^"), and the proposition is proved. 

Case 2°) Let us prove that the sum of (j30p is o(gi(n)^"'"*^'n^") by induction on the number 
I Q\ of classes of Q. Let us denote, for B' class of Q, Q(M ^ the partition of [2M] obtained 
from Q by linking classes B and {M}, and Q(M B^M + 1 B') the partition obtained 
from Q by linking classes B and {M}, and then the classes B' (which has become B U {M} if 
B' = B) and {M + 1} Let J(n) be the set of elements j of [n] such that D{N,n)j / 0. Note 
that, as noted to step II, | J(n)| < qi{n), and, by definition of D{N + l,n), J{n) is also the set 
of elements j of [n] such that D{N + 1, n)j ^ 0. 

• If Q has only one class, by boundedness hypothesis and since | < the sum of 
^ is 0{qi{n)). 

• Assume that the conclusion holds to ranks 1,...,|Q| — 1. 

- If neither {M} nor {M + 1} are singletons of Q, we apply the lemma, as for case 1°), with 
V = Q and TC = {h ^ [2M] ; k{h) = 1}. It leads to the conclusion in the same way as in case 
1°). 

- If exactly one of {M},{M + 1}, say {M}, is a singleton of Q, define V = Q - {{M}}, 
partition of the set [2M] — {M}. Note that the sum of ()30p . which is equal to the sum of (I3ip . 
is 

27V 

[A{M,n), n ^^nih,n) 

2N 2N 
27V 2Ar 

= ( m,n),)[ Y n^^.(^'^))-( E E 

3€J{n) i&[n]T h=l BeQ i(z\n]Q(M~'B h=l 

hf^M B^{M} ^ ' 

So by the induction hypothesis, for all B £ V,B ^ {^}) we have 

27V 

i6[„]Q(A/-.S) h=l 

So it suffices to prove that we have 

27V 

2+M'M" 



(32) Y MM,n)A( Y n =n^i(^)'^''''^ 



But by boundedness hypothesis, A{M, n)j = 0(| J(n)|) = 0{qi{n)). So it suffices to prove 

that the second sum in ([32]) is o{qi{n)^^'+^n'^'^"). Define H = {h e [2M] - {M} ; k{h) = 1}. By 
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definition of the A{h,n)'s, for all h £ TC, for all i S {qi{n) + 1, . . . ,n}, Ai(h,n) = 0. So by the 
lemma, the second sum in ()32p is 0(gi(n)l^^l^'P^/^nl^''^l^^w/2)^ Since, as noted above, for any 
partition with x singletons of a set S, one has \X\ — x/2 < \S\/2, we have 

irnl + \Vn\-^-^< (2M - 1)^ = M' + M" - Ifl, 
Moreover, since for all h € [2M] — {M}, k{h) = 1 implies h €z TC, 

\Vn\ - Pn/^ <\\{h^ [2M] - {M} ; Kh) / 1}| < M". 

Thus, since qi{n) < n, n^'Pnl-Pn/^ = ^FK|-PH/2-Af"^M" ^ 0(g^(n)l^«l-PH/2-Af"„M")_ jjence 
q^^n)\Pn\-Pn/2^\Pn\-Pn/2 = 0(g,(n)l^«l-¥+FH|-^-Af"^M") ^ 0{q,{n)'''-'^n''"), thus the 
conclusion holds. 

- If {M} and {M + 1} are singletons of Q, define P = Q - {{M}, {M + 1}}, partition of the 
set [2M] - {M, M + 1}. Note that the sum of ([30]), which is equal to the sum of ([31]), is 

2N 

^ [A{M,n)jA{M + l,n)j, ^ J] ^.J/i,n) 

i-i'6J(n) ie([n]-{j,j'})T' /i=l 

/i^A/+l 

Now, let us define Q' to be the partition of [2M] obtained from Q by linking classes {M} and 
{M + 1}. The sum of ()30p can be written 

2Af 2Ar 2Af 

A{M,n),A{M+l,n),>)[ U ^^dh,n))- n^*'^(^'™)-E E H^^'^^^'") 

ij'eJ{n) ie[n]P h=l i^[n]Q'h=l Tl ie[n]l^h=l 

where in the last sum, TZ runs over the set {Q(M B,M + 1 B') ; B,B' G Q}. So by the 
induction hypothesis, it suffices to prove that 

2N 

(33) ( AiM,n)jA{M + l,n),,)[Yl U A^ih,n)) = o {qr{nf+''n^") . 

j,j'eJ(n) iefnl'P h=l 

But by boundedness hypothesis, Ei,i'eJ(n) A{M,n)jA{M + l,n)j/ = 0(| J(n)|2) = 0(gi(n)2). So 

it suffices to prove that the second sum in ([33]) is o(gi(n)*^'n^^"). Define H := {h e [2M] - 
{M,M+l} ; k{h) = 1}. By definition of the ^(/i, n)'s, for all heH, for alH G {gi(n) + l, . . . ,n}, 
Ai{h,n) = 0. So by the lemma, the second sum in ()33p is 

Since, as noted above, for any partition X with x singletons of a set S, one has \X\ —x/2 < \S\/2, 
we have 

\Vn\ + \Vn\-^-^< (2M - 2)f2 = M' + M" - 1, 
Moreover, since for all h G [2M] - {M, M + 1}, k{h) = 1 implies heH, 

\Vn\-Pn/'^<l\{h^[m-{M,M + l};k{h)^l}\<M". 
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Thus, since qi{n) < n, 
Hence 

thus the conclusion holds. □ 
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